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Michael  Barnsley  made  two  visits  to  Wright  Patterson  to  discuss  opportu¬ 
nities  for  deployment  of  superfractals.  Based  on  collaborative  discussions  with 
Dr.  Ewing’s  team  and  Dan  Repperger,  and  feedback  from  the  presentation  made 
during  the  October  2008  trip,  he  made  a  second  presentation  in  March  2009,  at 
"Radar  Workshop  on  SuperResolution".  In  this  lecture  he  argued  that  it  is  fea¬ 
sible  to  use  superfractals  in  the  design  of  antennas.  This  would  have  advantages 
over  standard  fractal  designs:  (i)  implementation  of  E-variable  models,  whereby 
the  structure  of  the  antenna  is  of  several  characters  at  each  level;  (ii)  simplicity 
of  reconfiguration  to  different  designs,  for  testing.  The  approach  would  be  based 
on  the  use  of  fractal  homeomorphisms. 

A  number  of  papers  related  to  this  topic  were  written  during  the  period  of  the 
grant  and  are  listed  below.  The  published  versions  of  [1]  and  [2]  should  include 
explicit  acknowledgements  of  support  from  AOARD.  (Barnsley  also  published 
other  papers  during  this  period  but  these  were  not  related  to  the  topic.) 

A  high  level  programmer  with  Matlab  expertise  was  retained  to  investigate 
the  feasibility  of  making  an  image  recognition  system  based  on  fractal  homeo¬ 
morphisms,  in  a  Matlab  environment.  In  view  of  the  decision  to  focus  on  antenna 
design,  the  poor  performance  of  Matlab  code  in  comparision  with  Barnsley’s  ex¬ 
isting  code,  and  the  expense  of  continuing  to  work  with  Matlab,  this  direction 
was  discontinued. 

Some  progress  towards  an  intellectually  satisfying  model  for  living  systems, 
using  discrete  superfractals,  was  made. 
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THE  LIFE  AND  SURVIVAL  OF  MATHEMATICAL  IDEAS 


MICHAEL  F.  BARNSLEY 


Nature  and  evolution  provide  the  notion  of  a  creative  system:  a  core  stable  form 
(DNA),  a  fertile  environment,  a  determination  to  survive,  and  random  stimuli. 
Analogously,  the  mind  of  a  mathematician  provides  a  locus  for  creative  systems,  a 
place  where  mathematical  structures  live  and  evolve. 

According  to  the  King  James  version  of  Genesis,  "On  the  Fifth  day  ...  God 
created  great  whales."  But  Darwin  went  one  better;  at  the  end  of  his  masterwork, 
in  simple  beautiful  language,  he  proposes  that  what  was  created  was  a  creative 
system.  The  last  paragraph  of  Origin  of  Species  says: 

It  is  interesting  to  contemplate  a  tangled  bank,  clothed  with  many 
plants  of  many  kinds,  with  birds  singing  on  the  bushes,  with  various 
insects  flitting  about,  and  with  worms  crawling  through  the  damp 
earth,  and  to  reflect  that  these  elaborately  constructed  forms,  so 
different  from  each  other,  and  dependent  upon  each  other  in  so 
complex  a  manner,  have  all  been  produced  by  laws  acting  around 

us .  There  is  grandeur  in  this  view  of  life,  with  its  several 

powers,  having  been  originally  breathed  by  the  Creator  into  a  few 
forms  or  into  one;  and  that,  whilst  this  planet  has  gone  circling  on 
according  to  the  fixed  law  of  gravity,  from  so  simple  a  beginning 
endless  forms  most  beautiful  and  most  wonderful  have  been,  and 
are  being  evolved. 

Faced  with  the  extraordinary  richness  and  complexity  of  the  physical  observable 
universe,  of  which  we  are  part,  what  on  earth  can  a  mathematician  truly  create? 
The  answer  is:  a  vast  landscape  of  lovely  constructions,  born  for  the  first  time,  to 
live  on  in  the  realm  of  ideas.  For  the  realm  of  ideas  belongs  to  sentient  beings  such 
as  us:  whether  or  not  there  was  a  Creator,  it  is  certain  that  the  system  of  which 
we  are  part  is,  by  its  very  nature,  creative.  Our  genes  are  creative,  they  have  to  be, 
and  they  have  to  allow  creative  mutations.  They  must  be  stable  in  their  creativity. 
Their  creativity  is  the  well-spring  of  ours.  Not  only  must  our  genes,  through  the 
mechanisms  of  biology  and  random  mutation,  invent  new  viable  forms:  they  must 
be  prone  to  do  so. 

Our  mathematical  creativity  may  actually  be  initiated  by  random  events  at  the 
deepest  level,  after  deductive  reasoning,  consistencies,  experience,  and  even  intu¬ 
ition,  are  factored  out  of  the  process.  But  the  creative  mind  contains  something 
much  more  important  than  a  random  idea  generator;  it  provides  an  environment  in 
which  the  wild  seed  of  a  new  idea  is  given  a  chance  to  survive.  It  is  a  fertile  place. 
It  has  its  refugias  and  extinctions.  In  giving  credit  for  creativity  we  really  praise 
not  random  generation  but  the  determination  to  give  life  to  new  forms. 

But  when  does  a  newly  thought  up  mathematical  concept,  C,  survive?  Obviously, 
C  must  be  consistent  with  mathematics,  true,  correct,  etc.  But  I  believe  that  what 
causes  C  to  survive  in  the  minds  and  words  of  mathematicians  is  that  it  is,  itself. 
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a  creative  system.  For  now  I  will  resist  trying  to  give  a  precise  definition:  for  this 
young  notion  to  survive,  it  needs  to  be  adaptable.  Roughly,  I  mean  that  C  has 
the  following  attributes:  (I)  C  is  able  to  define  diverse  forms  and  structures,  call 
them  plants;  (II)  plants  possess  DNA;  (III)  C  is  stable  in  several  senses;  (IV)  C 
can  be  treated  from  diverse  mathematical  positions  (e.g.  topological,  geometrical, 
measure  theoretic,  algebraic) ;  (V)  C  is  highly  adaptable  and  can  be  translated  into 
the  languages  of  various  branches  of  science  and  engineering,  with  real  applications; 
(VI)  creative  systems  beget  new  creative  systems. 

For  more  that  two  thousand  years  the  key  forms  of  Euclid’s  Geometry  have 
survived,  shifted  in  importance,  and  evolved.  It  has  all  six  properties  of  a  creative 
system;  indeed  one  could  find  a  number  of  different  ways  of  defining  it  as  such.  Here 
is  the  one  that  I  like:  the  diverse  forms  and  structures  are  objects  such  as  lines, 
circles  and  other  constructions;  the  DNA  of  these  plants  are  formulas,  such  as  "the 
equation  for  a  straight  line",  provided  by  Descartes  analytic  geometry;  Euclidean 
geometry  can  be  treated  from  geometric,  algebraic  and  other  viewpoints;  the  topic  is 
stable  both  in  the  sense  that  nearby  DNA  yields  nearby  forms  and  in  the  sense  that 
small  changes  in  the  axioms  lead  to  new  viable  geometries;  it  has  adapted  to  many 
branches  of  science  and  engineering,  with  rich  applications;  and  Euclidean  geometry 
begat  projective  geometry  via  the  inclusion  of  the  line  at  infinity.  Alternatively  one 
might  describe  Euclid’s  geometry  more  abstractly  so  that  the  theorems  are  its  the 
diverse  structures  and  the  axioms  and  definitions  are  its  DNA. 

Dynamical  systems  [19]  and  eellular  automata  [35]  provide  two  recent  examples 
of  creative  systems.  I  mention  these  topics  because  each  has  an  obvious  visible 
public  aspect,  more  colorful  than  lines  and  circles  drawn  on  papyrus:  their  depth 
and  beauty  are  advertised  to  a  broad  audience  via  computer  graphics  represen¬ 
tations  of  some  of  their  flora.  They  are  alive  and  well,  not  only  in  the  minds  of 
mathematicians,  but  also  in  many  applications. 

In  your  own  mind  you  give  local  habitation  and  a  name  to  some  special  parts 
of  mathematics,  your  creative  system.  Since  this  note  is  a  personal  essay,  I  focus 
on  ideas  extracted  from  my  own  experience  and  research.  In  particular,  I  discuss 
iterated  function  systems  as  a  creative  system,  to  illustrate  connections  with  artistic 
creativity.  To  sharpen  the  presentation  I  focus  almost  exclusively  on  point-set 
topology  aspects.  While  the  specifics  of  iterated  function  systems  may  not  be 
familiar  to  you,  I  am  sure  that  the  mathematical  framework  is  similar  to  ones  that 
you  know. 

1.  Iterated  function  systems,  their  attractors,  and  their  DNA 

An  iterated  funetion  system  (IPS), 

T-=  (X;/i,...,/^), 

consists  of  a  complete  metric  space  X  together  with  a  finite  sequence  of  continuous 
functions,  {/„  :  X  ^  •  We  say  that  iF  is  a  contraetive  IFS  when  all  its  func¬ 

tions  are  contractions.  A  typical  IFS  ereative  system  may  consist  of  all  IFSs  whose 
functions  belong  to  a  restricted  family,  such  as  affine  or  bilinear  transformations 
acting  on 

Let  H  denote  the  set  of  nonempty  compact  subsets  of  X.  We  equip  El  with  the 
Hausdorff  metric,  so  that  it  is  a  complete  metric  space.  The  Hausdorff  distance 
between  two  points  in  HI  is  the  least  radius  such  that  either  set,  dilated  by  this 
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Figure  1.  Pictures  of  attractors  of  affine,  bilinear,  and  projective  IFSs. 


radius,  contains  the  other  set.  We  define  a  continuous  mapping  iF  :  HI  ^  H  by 

T{B)=UU{B) 

for  all  i?  e  H.  Note  that  we  use  the  same  symbol  T  for  the  IFS  and  for  the  mapping. 

Let  us  write  to  denote  the  composition  of  T  with  itself  k  times.  Then  we 
say  that  a  set  A  C  X  is  an  attractor  of  the  IFS  T  when  A  G  El  and  there  is  an  open 
neighborhood  M  oi  A  such  that 

lim  J^°^{B)  =  A 

k — >-oo 

for  all  i?  C  Af  with  i?  G  H.  Since  T  :  H  ^  HI  is  continuous  we  have  T  (A)  =  A. 
Notice  that  our  definition  of  attractor  is  topological:  in  the  language  of  dynamical 
systems,  A  is  a  strongly  stable  attractive  fixed-point  of  T . 

Our  first  theorem  provides  a  sufficient  condition  for  an  IFS  to  possess  an  attrac¬ 
tor,  one  of  the  "plants  of  many  kinds"  of  an  IFS  creative  system. 

Theorem  1.  [18]  Let  T  =  (X;  /i, ...,  /^r)  he  a  contractive  IFS.  Then  :  HI  ^  HI  is 
a  contraction,  and  hence,  by  Banach’s  contraction  theorem,  T  possesses  a  unique 
global  attractor. 

Attractors  of  IFSs  are  our  main  examples  of  the  diverse  forms  and  structures  of 
an  IFS  creative  system,  see  Figure  1,  as  in  attribute  (I).  An  affine  IFS  is  one  in 
which  the  mappings  are  affine  on  a  Euclidean  space.  Attractors  of  affine  IFSs  are 
the  bread-and-butter  sets  of  fractal  geometers,  such  as  Sierpinski  triangles,  twin- 
dragons,  Koch  curves.  Cantor  sets,  fractal  ferns,  and  so  on.  The  geometries  and 
topologies  of  these  attractors  are  so  rich,  fascinating  and  diverse  that  deep  papers 
get  written  about  a  single  species,  or  very  small  families  of  them! 

We  define  the  DNA  of  an  IFS  attractor  to  be  an  explicit  formula  for  the  IFS.  We 
refer  to  the  DNA  of  an  IFS  attractor  as  an  IFS  eode.  The  DNA  for  the  canonical 
Cantor  set  is  (K;  fi{x)  =  a:/3,  f2{x)  =  (a:  -I-  2)  /3);  these  few  symbols  and  their  con¬ 
text  define  a  non-denumerable  set  of  Lebesque  measure  zero.  Similarly,  the  DNA 
for  the  Sierpinski  triangle  is  (K^;  /i(x,  y)  =  (x/2,  y/2),  /2(a:,  y)  =  {x/2  +  1/2,  y/2), 
fz{x,y)  =  {xl2,yl2  +  1/2)).  Here  a  curve  whose  points  are  all  branch  points  is 
captured  in  a  short  strand  of  symbols.  Other  simple  IFS  codes  provide  DNA  for 
classical  objects  such  as  arcs  of  parabolas,  line  segments,  triangles,  and  circles. 

How  do  the  individual  numbers  in  IFS  codes  relate  to  the  properties  of  the 
attractors  that  they  define?  Similarly,  we  might  ask  about  the  relationship  between 
the  DNA  of  a  biological  plant  and  the  plant  itself,  the  details  of  its  leaf  shapes,  the 
structure  of  its  vascular  bundles,  and  so  on. 
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2.  When  does  an  aefine  IFS  possess  an  attractor? 


In  discussing  this  seemingly  simple  question  we  reveal  how  IFS  theory  is  subtle 
and  leads  into  applications,  as  in  attribute  (V).  We  characterize  both  geometrically 
and  metrically,  see  attribute  (IV),  those  affine  IFSs  that  possess  attribute  (I). 

Intuition  incorrectly  suggests  that  the  answer  to  our  question  is:  if  the  magni¬ 
tudes  of  all  of  the  eigenvalues  of  the  linear  parts  of  the  maps  of  the  IFS  are  less 
than  one,  then  the  affine  IFS  has  an  attractor.  The  situation  seems  to  be  analogous 
to  the  situation  for  discrete  dynamical  systems,  [17]  Proposition,  p.279,  where  an 
affine  map  has  an  attractive  fixed  point  if  and  only  if  the  norm  of  the  linear  part  is 
less  than  one.  But  the  situation  is  not  analogous.  Consider  for  example  the  IFS 

=  {2y,-x/3),f2ix,y)  =  {-y/3,2x)). 

The  point  O  =  (0,0)  is  an  attractive  fixed  point  for  both  fi  and  /2,  because 
ff'^{x,y)  =  /2°^"(x,y)  =  (-2/3)^  {x,y) ,  ,  y)  =  (-2/3)"/i(x,  y),  and 

y.o(2n-i-i)j-^^ =  (— 2/3)"/2(a;, y)  for  all  points  {x,y)  in  But  {O}  is  not  an 
attractor  for  the  IFS  because  (/i  o  /2)”(a;,y)  =  (4"ai,y/9”)  implies  that  O  is  an 
unstable  fixed  point. 

The  following  theorem  contains  an  answer  to  our  question,  and  an  affine  IFS 
version  of  the  converse  to  Banach’s  contraction  theorem  [13].  For  us,  most  impor¬ 
tantly,  it  provides  both  a  metric  and  a  geometrical  characterization  of  viable  affine 
IFS  codes.  See  also  Berger  and  Wang,  [9],  and  Daubechies  and  Lagarias,  [10]. 

Theorem  2.  [1]  If  T  =  (M-^; /i, ..., /at)  is  an  affine  IFS  then  the  following  state¬ 
ments  are  equivalent. 

(1)  T  possesses  an  attractor. 

(2)  There  is  a  metric,  Lipshitz  equivalent  to  the  Euclidean  metric,  with  respect 
to  which  each  fn  is  a  contraction. 

(3)  There  is  a  closed  bounded  set  K  C  ,  whose  affine  hull  is  such  that 
T  is  non-antipodal  with  respect  to  K . 


Briefly,  let  me  explain  the  terminology.  We  say  that  two  metrics  di  (•,•)  and 
(i2  (•,  •)  on  R^  are  Lipshitz  equivalent  when  there  is  a  constant  C  >  1  such  that 
di{x,y)/C  <  d2{x,y)  <  Cdi{x,y)  for  all  x,y  G  Given  any  closed  bounded 

set  K  in  R^ ,  whose  affine  hull  is  R^ ,  and  any  u  €  S^~^,  the  unit  sphere  in 
R^ ,  let  {TCujTC-u}  be  the  unique  pair  of  distinct  support  hyperplanes  of  K  per¬ 
pendicular  to  u,  see  [27]  p.l4.  Then  the  set  of  antipodal  pairs  of  points  of  K  is 
K'  :=  {{a,  a'}  :  a  €  TLu  H  dK,  a'  G  TL-u  C  dK,  u  G  S^~^}  where  dK  denotes  the 
boundary  of  K.  We  say  that  an  IFS  T  is  non-antipodal  with  respect  to  K  when  each 
of  its  functions  takes  K  into  itself  but  maps  no  antipodal  pair  of  points  of  K  to  an 
antipodal  pair  of  points  of  K.  We  denote  the  latter  condition  by  T  [K')  n  iG'  =  0. 

Part  of  the  proof  of  Theorem  2  relies  on  the  observation  that  if  /C  C  R^  is  a 
convex  body  (think  of  K.  as  the  convex  hull  of  K  in  (3))  then  we  can  define  a  metric 
dic{-,-)  on  R^ ,  Lipshitz  equivalent  to  the  Euclidean  metric,  by 


dK{x,y)  =  inf 


\\x-y\\ 

\\l-m\\ 


I,  toG/C,  I  —  m  =  a(x 


y),a& 


for  all  X  y,  where  \\x  —  y\\  denotes  the  Euclidean  distance  from  x  to  y  in  R^ . 
One  shows  that,  if  an  affine  map  /„  is  non-antipodal  with  respect  to  fC,  then  it  is  a 
contraction  with  respect  to  djc.  In  fact  dfc  is,  up  to  a  constant  factor,  a  Minkowski 
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Figure  2.  The  left-hand  panel  shows  a  black  fern  image  within  a 
grimy  window;  overlayed  upon  it  are  four  affine  transformations  of 
the  window  and  the  fern,  with  the  transformed  fern  images  shown 
in  green.  The  goal  has  been  to  approximate  the  original  (black) 
fern  with  affinely  transformed  (green)  copies  of  itself.  The  rectan¬ 
gular  window  has  been  mapped  non-antipodally  upon  itself.  So, 
by  Theorem  2,  there  exists  a  metric  such  that  the  associated  affine 
IFS  T  is  contractive.  The  original  fern  (black)  and  the  attractor 
(red)  of  T  are  shown  in  the  right-hand  panel. 


metric  [30]  associated  with  the  symmetric  convex  body  defined  by  the  Minkowski 
difference  1C  —  1C. 

Theorem  2  provides  a  means  for  defining  viable  IFS  codes  (DNA)  and  is  useful 
in  the  design  of  a  two-dimensional  affine  IFS  whose  attractor  approximates  a  given 
target  set  T  C  Typical  IFS  software  for  this  purpose  exhibits  a  convex  window 
/C,  containing  a  picture  of  T,  on  a  digital  computer  display.  A  set  of  affine  maps 
is  introduced,  thereby  defining  an  IFS  T .  The  maps  are  adjusted  using  interactive 
pictures  of  T  (/C)  and  T  (T).  If  we  ensure  that  T  {JC)  C  JC  and  T  {K.')  fl  /C'  =  0, 
then  T  possesses  a  unique  attractor  A  C  1C.  An  example  is  illustrated  in  Figure  2. 

Now  we  are  in  a  position  to  explain  a  stability  relationship  between  IFS  codes 
and  attractors,  and  thus  to  exhibit  a  form  of  stability  as  in  attribute  (III).  We 
can  control  the  (Hausdorff)  distance  hjc,  which  depends  on  djc,  between  A  and  T 
because  it  depends  continuously  on  the  distance  between  fF  (T)  and  T.  Indeed,  the 
collage  theorem  [2]  states  that 

/.UAdsMAAT) 

where  0  <  A  <  I  is  a  Lipshitz  constant  for  fF  :  H  ^  H,  for  example  the  maximum 
of  a  set  of  contractivity  factors  of  the  /„s  with  respect  to  die.  Notice  that  this  rela¬ 
tionship  says  nothing  about  the  topological  structure  of  an  attractor:  it  comments 
only  on  its  approximate  shape. 

The  collage  theorem  expresses  one  kind  of  stability  for  the  IFS  creative  system, 
as  in  attribute  (III):  small  changes  in  the  IFS  code  of  a  contractive  IFS  lead  to  small 
changes  in  the  shape  of  the  attractor.  Indeed,  this  realization  played  a  role  in  the 
development  of  fractal  image  compression  [2].  In  this  development  several  things 
occurred.  First,  affine  IFS  theory  adapted  to  a  digital  environment,  illustrating  a 
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component  of  attribute  (V).  Since  contractive  affine  IFSs  do  not  in  general  translate 
to  contractive  discrete  operators,  [28],  new  theory  had  to  be  developed.  (See  for 
example  [16].)  This  illustrates  stability  of  a  second  kind,  as  required  in  attribute 
(III):  the  underlying  ideas  are  robust  relative  to  structural  changes  in  the  creative 
system.  Finally,  a  real  application  was  the  result,  as  required  by  attribute  (V). 


3.  Projective  and  Bilinear  IFSs 


We  will  also  use  both  projective  and  bilinear  IFSs  for  creative  applications  that 
control  the  shape  and  topology  of  attractors,  and  transformations  between  attrac¬ 
tors.  Both  are  generalizations  of  two-dimensional  affine  IFSs.  Both  can  be  ex¬ 
pressed  with  relatively  succinct  IFS  codes,  yet  have  more  degrees  of  freedom  than 
affines.  Another  such  family  of  IFSs  is  provided  by  the  Mobius  transformations  on 
C  U  {oo}.  The  availability  of  a  rich  selection  of  accessible  examples  is  a  valuable 
attribute  for  the  survival  of  a  mathematical  idea. 

In  two  dimensions,  the  functions  of  a  projective  IFS  are  represented  in  the  form 


(3.1) 


fn{x,y)  =  ( 


anX  -I-  bny  +  Cn 
gnX  +  h„y  +  jn 


dnX  +  Cny  +  kn 

QnX  +  hny  +  jn 


), 


where  the  coefficients  are  real  numbers.  A  similar  result  to  Theorem  2  applies 
to  such  projective  transformations  restricted  to  a  judiciously  chosen  convex  body, 
with  the  associated  Hilbert  metric,  [12]  p.l05,  used  in  place  of  the  generalized 
Minkowski  metric  mentioned  above.  Specihcally,  let  T  denote  a  projective  IFS  of 
the  form  (/C°;  /i, ...,  /at),  where  /C°  is  the  interior  of  a  convex  body  /C  C  such 
that  T  (/C)  C  K.° .  The  associated  Hilbert  metric  dn  is  defined  on  /C°  by 


dnix,  y)  =  In  \TZ  (x,  y;  a,  6)|  for  all  x,y  €  IC°  with  x  ^  y, 


where  TZ  {x,  y;  a,  b)  =  (|6  —  x\l  \x  —  a|)  /  (|6  —  y\/  \y  —  a|)  denotes  the  cross  ratio  be¬ 
tween  X,  y  and  the  two  intersection  points  a,  b  of  the  straight  line  through  x,  y  with 
the  boundary  /C.  You  might  like  to  verify  that  iF  is  a  contractive  IFS  with  re¬ 
spect  to  dn,  using  the  fact  that  projective  transformations  preserve  cross  ratios. 
So  projective  IFSs  can  be  used  in  applications  in  nearly  the  same  way  as  affine 
systems. 

To  describe  bilinear  transformations,  let  TZ  =  [0, 1]^  C  denote  the  unit  square, 
with  vertices  A  =  (0, 0),  B  =  (1, 0),  C  =  (1, 1),  D  =  (0, 1).  Let  P,  Q,  R,  S  denote,  in 
cyclic  order,  the  successive  vertices  of  a  possibly  degenerate  quadrilateral.  Then  we 
uniquely  dehne  a  bilinear  function  B  :TZ  ^  TZ  such  that  B{ABCD)  =  PQRS  by 


(3.2)  B{x,  y)  =  P  +  x{Q  -  P)  +  y{S  -  P)  +  xy{R  P  P  -  Q  -  S). 


This  transformation  acts  affinely  on  any  straight  line  that  is  parallel  to  either  the 
x-axis  or  the  y-axis.  For  example,  if  B\ab  '■  AB  PQ  is  the  restriction  of  B  to 
AB,  and  if  Q  :  ^  is  the  affine  function  dehned  by  Q{x,  y)  =  P  +  x{Q  —  P)  + 

y{S  —  P),  then  Q\ab  =  B\ab-  Because  of  this  "affine  on  the  boundary"  property, 
bilinear  functions  are  well  suited  to  the  construction  of  fractal  homeomorphisms, 
as  we  will  see.  Sufficient  conditions  under  which  there  exists  a  metric  with  respect 
to  which  a  given  bilinear  transformation  is  contractive  are  given  in  [6].  A  bilinear 
IFS  has  an  attractor  when  its  IFS  code  is  close  enough  (in  an  appropriate  metric) 
to  the  IFS  code  of  an  affine  IFS  that  has  an  attractor. 

An  example  of  a  geometrical  configuration  of  quadrilaterals  that  gives  rise  to 
both  a  projective  and  a  bilinear  IFS  is  illustrated  in  Figure  3.  In  either  case  we 
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Figure  3.  The  four  quadrilaterals  IE  AH,  lEBF,  IGCE,  IGDH, 
define  both  an  projective  E  and  a  bilinear  IFS  Q  both  of  which  have 
a  unique  attractor,  the  filled  rectangle  with  vertices  at  ABGD; 
but  the  address  structures  are  different.  In  the  left  panel  and 
right  panels,  the  attractors  of  T  and  Q  respectively  have  been 
rendered,  so  that  points  with  the  same  address  are  the  same  shade 
of  grey.  The  address  structure  is  independent  of  E,  E,  G,  H,  I,  in 
the  bilinear  case  but  not  in  the  projective  case.  (Hint:  compare 
how  lines  meet  the  line  IE.) 


define  the  IFS  to  be  ilZ]  /i,  /2,  /s,  /t)  where 

fi{ABGD)  =  IE  AH,  (ABGD)  =  lEBF,  /a  (ABGD)  =  IGGF,  fsiABGD)  =  IGDH, 

where  the  first  expression  means  /i  (H)  =  I,  fi  (B)  =  E,  fi  (G)  =  A,  fi  (D)  =  H. 

With  few  constraints  each  IFS  is  contractive  with  respect  to  a  metric  that  is  Lipshitz 
equivalent  to  the  Euclidean  metric,  with  attractor  equal  to  the  hlled  rectangle 
ABGD.  But  there  is  an  important  difference:  the  bilinear  family  provides  a  family 
of  homeomorphisms  on  IZ,  with  applications  to  photographic  art,  attribute  (V), 
while  the  projective  family  does  not,  as  we  will  see. 

4.  The  Chaos  Game 

How  do  we  compute  approximate  attractors  in  a  digital  environment?  Algo¬ 
rithms  based  on  direct  discretization  of  the  expression  A  =  limfc^oo  have 

high  memory  requirements  and  tend  to  be  inaccurate,  [28].  The  availability  of  a 
simple  algorithm  that  is  fast  and  accurate,  for  the  types  of  IFS  that  we  discuss,  has 
played  an  important  role  in  the  survival  of  the  IFS  creative  system.  The  following 
algorithm,  known  as  the  chaos  game,  was  described  to  a  wide  audience  in  Byte 
Magazine  in  1988,  and  successfully  dispersed  IFS  codes  to  the  computer  science 
community.  It  helped  ensure  that  the  IFS  creative  system  would  have  attribute 
(V). 

Define  a  random  orbit  {xk}^^Q  of  a  point  Xq  G  X  under  T  =  (X; /i, ..., /jv) 
by  Xk  =  /crj, (xfc-i),  where  Ok  G  {1,2,  ...,A^}  is  chosen,  independently  of  all  other 
choices,  by  rolling  an  iV-sided  die.  The  random  orbit  is  associated  with  cr  G  H.  If 
the  underlying  space  is  two-dimensional  and  F  is  contractive,  then  it  is  probable 
that  a  picture  of  the  attractor  of  T ,  accurate  to  within  viewing  resolution,  will  be 
obtained  by  plotting  on  a  digital  display  device. 

Why  does  this  Markov  chain  Monte  Carlo  algorithm  work?  The  following  theo¬ 
rem,  implicit  in  [8],  tells  us  how  we  can  think  of  the  attractor  of  a  contractive  IFS 
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(I)  Initialize  z«0.  y«0. 

(ii)  For  n«>l  10  2500.  do  step*  (ni)-(vn). 

(iii)  Okxmc  a  to  be  one  of  the  numben  1 . 

2.  ..  m.  with  probability  p.. 

(iv)  A|^y  the  (rantfomuilion  to  the 
point  U.yiioobuin  (i.y). 

(V)  Set  (x.y)  equal  to  the  nr*  point  x«i, 

ymy 

(vi) lfn>IO.  ploKc.yl 

(vii)  Loop 


Figure  4.  The  original  pseudocode  from  Byte  Magazine  (Jan 
1988)  for  implementing  Chaos  Game  algorithm  to  obtain  an  image 
of  the  attractor  of  an  IFS  on  Notice  the  small  number  of  it¬ 
erations  used!  Nowdays,  usually,  I  use  10^  iterations,  and  discard 
the  first  thousand  points.  On  the  right  is  a  sketch  of  a  2-variable 
tree  obtained  by  a  generalization  of  the  chaos  game. 


Figure  5.  From  left  to  right  this  picture  shows:  the  result  of 
9000  iterations  of  the  chaos  game  algorithm  applied  to  a  projective 
IFS;  the  result  of  10^  iterations;  a  small  picture  of  a  flower;  and  a 
rendered  close-up  of  the  attractor.  In  the  latter  image  the  colours 
were  obtained  with  the  aid  of  a  fractal  transformation  from  the 
attractor  to  the  small  picture  of  the  yellow  flower. 


as  being  the  w-limit  set  of  almost  any  random  orbit.  A  direct  proof  can  be  found 
in  [31]. 

Theorem  3.  Let  be  a  random  orbit  of  a  contractive  IFS.  With  probability 

one 

lim  {xk}  =  A, 

iC  — >-oo 

where  the  limit  is  taken  with  respect  to  the  Hausdorff  metric. 

Pictures,  calculated  using  the  chaos  game  algorithm,  of  the  attractor  of  a  pro¬ 
jective  IFS  (7?.;  /i,  /2,  /s,  fi)  are  shown  in  the  left-most  two  panels  of  Figure  5.  You 
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can  visualize  an  approximation  to  the  stationary  probability  measure  of  the  sto¬ 
chastic  process,  implicit  in  the  chaos  game,  in  the  left-hand  image.  This  measure 
depends  on  the  strictly  positive  probabilities  associated  with  the  maps,  but  its  sup¬ 
port,  the  attractor,  does  not!  Measure  theory  aspects  of  IFS  are  not  considered  in 
this  article,  but  it  is  nice  to  see  one  way  in  which  the  topic  arises. 


5.  Addresses,  and  transformations  between  attractors 

In  this  section  we  deepen  our  understanding  of  attractors.  We  discover  informa¬ 
tion  about  the  relationship  between  IFS  codes  and  the  topology  of  attractors,  and 
the  relationships  between  different  attractors.  This  information  helps  to  classify 
our  diverse  plants,  attribute  (I),  and  leads  into  real  applications  in  art  and  biology, 
attribute  (V).  Good  book-keeping  is  the  key. 

The  space  =  {1,  2, ...,  N}°°  with  the  product  topology  plays  a  fundamental  role 
in  IFS  theory  and  in  this  article.  We  write  a  =  aia2---  to  denote  a  typical  element 
of  n.  We  will  use  the  notation  '■=  Ui  o  °  ■■■  °  Uk,  ^\k  =  CTi(J2...(Tfe, 

and  /crio-2...o'fc  =  fa\k  for  any  a  efl  and  fc  =  I,  2, ...  . 

The  following  theorem  suggests  our  plants  can  have  intricate  topological  struc¬ 
tures  and  suggests  that  symbolic  dynamics  are  involved,  thereby  adding  a  lusher 
interpretation  of  attribute  (I). 

Theorem  4.  [18]  Let  T  =  (X; /i, ..., /^v)  he  a  contractive  IFS,  with  attractor  A. 
Let  X  G  X.  A  continuous  surjection  tt  :  LI  ^  A,  independent  of  x,  is  well-defined 
by  TT  (a)  —  limfc^oo  fa\k{x);  the  convergence  is  uniform  for  (cr,x)  G  Lt  x  B,  for  any 
B  gU. 

The  set  of  addresses  of  a  point  x  G  A  is  defined  to  be  the  set  7r“^  (x),  and 
defines  an  equivalence  relation  ~  on  Ll.  For  example,  the  attractor  of  the  IFS 
(M;  /i(x)  =  x/2,  /2(x)  =  x/2  -I-  1/2)  is  the  closed  interval  [0,  Ij.  You  may  check  that 
7r-i(0)  =  {T  :=  nil...}  ,7r-i(l)  =  {2},  7r-i(l/2)  =  {12, 2T},  and  7r-i(l/3)  = 
{12}.  Some  points  of  an  attractor  have  one  address  while  others  have  multiple 
distinct  addresses.  The  topology  on  A  is  the  identification  topology  on  Ll  induced 
by  the  continuous  map  :  LI  ^  A.  In  this  paper  we  refer  to  the  set  of  equivalence 
classess  induced  by  ~  on  17  as  the  address  structure  of  the  IFS.  Figure  3  contrasts 
the  address  structures  of  a  corresponding  pair  of  bilinear  and  projective  IFSs. 

We  can  think  of  the  topology  of  an  attractor  A  as  being  that  of  17  with  all  points 
in  each  equivalence  class  glued  together;  that  is,  A  is  homeomorphic  to  17/ 
Simple  examples  demonstrate  that  the  address  structure  can  change  in  complicated 
ways  when  a  single  parameter  is  varied:  the  topologies  of  attractors,  our  plants,  in 
contrast  to  their  shapes,  do  not  in  general  depend  continuously  on  their  IFS  codes. 
By  restricting  to  appropriate  families  of  projective  or  bilinear  IFSs,  with  known 
address  structures,  control  of  the  topology  of  attractors  becomes  feasible. 

A  point  on  an  attractor  may  have  multiple  addresses.  We  select  the  "top"  address 
to  provide  a  unique  assignment;  the  top  address  is  the  one  closest  to  I  =  fill...  in 
lexographic  ordering.  Each  element  of  the  address  structure  of  an  IFS  is  represented 
by  a  unique  point  in  17.  This  choice  is  serendipitous,  because  the  resulting  set  of 
addresses,  called  the  tops  space,  is  shift  invariant,  and  so  yields  a  link  between 
our  plants,  symbolic  dynamics,  and  information  theory,  attribute  (IV),  see  [5]  and 
references  therein. 
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Figure  6.  Three  renderings  of  a  close-up  of  the  attractor  in  Fig¬ 
ure  5  computed  using  a  coupled  version  of  the  chaos  game  and 
a  fractal  transformation.  On  the  left  the  computation  has  been 
stopped  early,  yielding  a  misty  effect.  Different  aspects  of  the  frac¬ 
tal  transformation  are  revealed  by  applying  it  to  different  pictures. 


We  define  a  natural  map  from  an  attractor  Aj:  of  an  IFS  T  to  the  attractor  Aq 
of  an  IFS  Q,  each  with  the  same  number  of  maps,  by  assigning  to  each  point  of  Ajr 
the  point  of  Ag  whose  set  of  addresses  includes  the  top  address  of  the  point  in  Ajr. 
This  provides  a  map  Tj^g  :  Ajr  Ag  called  a  fractal  transformation.  When  the 
address  structures  of  Ajr  and  Ag  are  the  same,  this  map  is  a  homeomorphism.  Since 
fractal  transformations  can  be  readily  computed  by  means  of  a  coupled  version  of 
the  chaos  game,  applications  to  art  and  geometric  modelling  become  feasible  and 
the  IFS  creative  system  tests  new  forms  and  environments,  attribute  (V). 

Let  7?.  C  denote  a  filled  unit  square.  Let  p  :  TZ  ^  C  he  d.  picture  (function); 
that  is,  p  is  a  mapping  from  TZ  into  a  color  space  C.  A  color  space  is  a  set  of  points 
each  of  which  is  associated  with  a  unique  color.  In  computer  graphics  a  typical 
color  space  is  C  =  {0, 1, ...,  255}^,  where  the  coordinates  of  a  point  represent  digital 
values  of  red,  green  and  blue.  The  graph  of  a  picture  function  may  be  represented 
by  a  colorful  picture  supported  on  TZ.  Next  time  you  see  a  picture  hanging  on  a 
wall,  imagine  that  it  is  instead,  an  abstraction,  a  graph  of  a  picture  function.  More 
generally,  we  allow  the  domain  of  a  picture  function  to  be  an  arbitrary  subset  of 
M2. 

For  example,  in  the  right-hand  image  in  Figure  5  we  have  rendered  the  graph 
of  p  :  Ayr  C  obtained  by  choosing  p  to  correspond  to  the  picture  of  the  yellow 
flower,  Q  to  be  an  affine  IFS  such  that  {pn  {TZ)}'^^i  is  a  set  of  rectangular  tiles 
with  Lipn  (TZ)  =  TZ,  and  T  to  be  the  projective  IFS  whose  attractor  is  illustrated  in 
black.  In  Figure  6  we  show  other  renderings  of  a  portion  of  the  attractor,  obtained 
by  changing  the  picture  p  and,  in  the  left-hand  image,  by  stopping  the  chaos  game 
algorithm  early-the  pictures  are  all  computed  by  using  a  coupled  variant  of  the 
chaos  game.  Since  p  =  po  Tjrg  one  can  infer  something  about  the  nature  of  fractal 
transformations  looking  at  such  pictures.  By  panning  the  source  picture  p  it  is 
possible  to  make  fascinating  video  sequences  of  images.  You  can  see  some  yourself 
with  the  aid  of  SF  Video  Shop  [32].  In  the  present  example  you  would  quickly  infer 
that  Tjrg  is  not  continuous  but  that  it  is  not  far  from  being  so:  it  may  be  continuous 
except  across  a  countable  set  of  arcs. 
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Figure  7.  (i)  The  adjustable  points  a,  b,  c  used  to  define  a  family 
of  afRne  iterated  function  systems  =  (M^;  /i,  /2,  /s,  /t)  with  con¬ 
stant  address  structure;  (ii)  a  picture  supported  on  the  attractor  of 
an  IFS  belonging  to  the  family;  (iii)  the  same  picture,  transformed 
under  a  fractal  homeomorphism  of  the  form  Tjrg  =  7rg  o  tjf  where 
T  and  Q  are  two  IFSs  belonging  to  the  family. 


6.  The  art  of  fractal  homeomorphism 

Affine  and  bilinear  iterated  function  systems  can  be  used  to  provide  a  wide 
variety  of  parameterized  families  of  homeomorphisms  on  two-dimensional  regions 
with  polygonal  boundaries  such  as  triangles  and  quadrilaterals.  We  use  them  to 
illustrate  the  application  of  the  IFS  creative  system  to  a  new  artform.  In  effect  this 
application  is  itself  a  new  creative  system,  for  artists.  This  provides  an  illustration 
of  attribute  (VI):  creative  systems  beget  creative  systems. 

For  example,  let  A,  B,  and  C  denote  three  non-colinear  points  in  K^.  Let  c  denote 
a  point  on  the  line  segment  AB,  let  a  denote  a  point  on  the  line  segment  BC,  and 
let  b  denote  a  point  on  the  line  segment  CA,  such  that  {a,b,c}  fl  {A,B,C}  =  0; 
see  Figure  7(i).  Let  T  =  (M^;  /i,  /2,  /a,  fi)  be  the  unique  afhne  IFS  such  that 

fi{ABC)  =  caB,  f2{ABC)  =  Cab,  f^{ABC)  =  cAb,  and  /^{ABC)  =  cab, 

where  we  mean  for  example  that  fi  maps  A  to  c,  B  to  a,  and  C  to  B,  see  Figure 
7(i).  For  reference,  let  us  write  T  =  where  a  =  \Bc\/\AB\,j5  =  \Ca\/\BC\, 

and  7  =  lA^I/ICAj.  The  attractor  of  is  the  filled  triangle  T  with  vertices 

at  A,  B,  and  C.  Then  tFa,p,-y  is  contractive  IFS,  for  each  ia,l3,j)  e  (0, 1)3,  with 
respect  to  a  metric  that  is  Lipshitz  equivalent  to  the  Euclidean  metric,  by  Theorem 
2.  Its  address  structure  Ca,i5,'y  is  independent  of  a,  13,^,  see  [5]  section  8.1. 

Figure  7(ii)  illustrates  a  picture  p  :  T  ^  C;  it  depicts  fallen  autumn  leaves.  Fig¬ 
ure  7(iii)  illustrates  the  picture  p  =  poTj^g,  namely  the  result  of  applying  the  home¬ 
omorphism  Tyrg  to  the  picture  p,  where  T  =  ^^0  45  0,45.0.45  and  Q  =  .7^0.55.0.55,0.55  ■ 
The  transformation  in  this  example  is  area-preserving  because  corresponding  tiles 
have  equal  areas. 

A  similar  result  applies  to  families  of  bilinear  IFSs.  For  example  Figure  3  de¬ 
fines  a  family  of  bilinear  IFSs,  T^,  parameterized  by  the  vector  of  points  v  = 
{E,  F,G,  H,  I).  This  family  has  constant  address  structure  for  all  values  of  v  for 
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Figure  8.  Example  of  a  fractal  homeomorphism  generated  by  two 
IFSs  of  four  bilinear  transformations.  The  attractor  of  each  trans¬ 
formation  is  the  unit  square. 


Figure  9.  Before  (left)  and  after  a  fractal  homeomorphism. 


which  T-u  is  contractive  and  can  thus  be  used  to  provide  a  family  of  homeomor- 
phisms  'Tv,w  :  7?.  ^  7?..  An  illustration  of  the  action  of  'Tv,w  on  a  picture  of  Aus¬ 
tralian  heather  is  given  in  Figure  8.  In  this  case  the  parameters  v  and  w  both 
correspond  to  affine  IFSs.  What  is  remarkable  in  this  case,  and  many  like  it,  is  that 
the  transformed  picture  looks  so  realistic.  Can  you  tell  which  is  the  original? 

Figure  9  illustrates  a  homeomorphic  fractal  transformation  generated  by  a  pair  of 
bilinear  IFSs  on  7?..  In  this  case  iV  =  12.  The  original  image  is  a  digital  photograph 
of  a  lemon  tree  and  wallflowers  in  my  garden  in  Canberra.  The  final  image  was 
printed  out  on  thick  acid-free  rag  paper  by  a  professional  printing  company,  using 
vivid  pigment  inks,  at  a  width  of  approximately  5ft  and  a  height  of  3ft  Gins.  It 
represents  a  fusion  of  the  colors  of  nature  and  mathematics;  it  provokes  wonder  in 
me,  a  sense  of  the  pristine  and  inviolate,  a  yearning  to  look  and  look  ever  closer, 
see  Figure  10. 

I  have  used  such  extraordinary  transformations  to  generate  works  for  three  suc¬ 
cessful  (most  of  the  pictures  are  sold)  art  shows,  in  Canberra  (Australia,  April 
2008),  in  Bellingham  (Washington  State,  April  2008),  and  in  Gainesville  (Florida, 
March  2009). 
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Figure  10.  Detail  from  Figure  9  showing  not  only  the  vibrant 
colors  of  nature,  but  also  the  wide  range  of  stretching  and  squeezing 
achieved  in  this  relatively  simple  fractal  transformation. 


7.  SUPERFRACTALS 

In  this  section  we  illustrate  attribute  (VI) .  We  show  how  IFS  theory  begets  a  new 
creative  system  via  a  higher  level  of  abstraction.  The  new  framework  is  suitable 
for  mathematical  modelling  of  the  geometry  of  a  multitude  of  naturally  occurring, 
readily  observable  structures.  It  also  has  applications  to  the  visual  arts. 

The  new  system  has  some  remarkable  properties.  Its  attractor  is  a  set  of  in¬ 
terrelated  sets  that  can  be  sampled  by  a  variant  of  the  chaos  game  algorithm,  as 
illustrated  in  Figure  12.  This  algorithm  is  born  fully  formed  and  is  the  key  to 
applications.  The  geometry  and  topology  of  the  interrelated  sets  can  be  controlled 
when  appropriate  generalized  IFSs  are  used.  In  particular,  through  the  concepts 
of  V-variability  [7]  and  superfractals,  we  are  able  to  form  a  practical  bridge  be¬ 
tween  deterministic  fractals  (such  as  some  of  the  IFS  attractors  in  previous  sections) 
and  random  fractal  objects  (such  as  statistically  self  similar  curves  that  represent 
Brownian  motion). 

7.1.  V— variability.  Here  is  a  biological  way  to  think  of  "V- variability".  Imagine 
a  tree  that  grows  with  this  property.  If  you  were  to  break  off  all  of  the  branches  of 
any  one  generation  and  classify  them,  you  would  find  that  they  were  of  at  most  V 
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Figure  11.  An  element  of  a  2-variable  superfractal  is  shown  on 
the  left.  If  you  look  closely  at  it,  you  will  see  that  it  is  made  of 
exactly  two  distinct(up  to  translation,  reflection  and  rotation)  sub¬ 
objects  of  half  the  linear  dimension.  (The  two  objects  immediately 
to  the  right  of  the  first  one.)  And,  if  you  look  even  closer,  you  will 
see  that  it  is  made  of  two  sub-objects  of  one-quarter  the  linear 
dimension,  and  so  on. 


different  types.  By  "generation"  I  mean  that  you  are  able  to  think  of  the  tree  as 
having  older  and  younger  branches,  that  is,  some  that  started  to  grow  during  year 
one,  sub-branches  that  began  during  year  two  and  so  on.  The  tree  is  very  old.  By 
"type"  I  mean  something  like  "belongs  to  a  particular  conjugacy  class".  The  type 
may  change  from  generation  to  generation,  but  the  number  V  is  fixed  and  as  small 
as  possible.  Then  I  will  call  the  imagined  tree  "F-variable".  Figure  4  includes 
an  illustration  of  a  2-variable  tree,  where  the  younger  branches  start  higher  up 
the  tree.  Again,  consider  a  population  of  annual  plants  belonging  to  a  species  that 
admits  S  distinct  possible  genotypes.  If  the  number  of  distinct  genotypes  in  each 
generation  is  bounded  above  by  V  then  (in  circumstances  where  V  is  significantly 
smaller  than  S)  I  would  call  this  population  "F-variable".  But  the  mathematical 
definition  relates  to  a  property  of  attractors  of  certain  IFSs.  Figure  11  illustrates 
a  2-variable  fractal  subset  of  the  Euclidean  plane:  it  is  a  union  of  two  tiles  of  half 
its  size:  it  is  also  a  union  of  at  most  two  tiles  of  a  quarter  its  size,  and  so  on. 

Let  T  =  (X; /i, ..., /tv)  be  an  IFS  of  functions  /„  that  are  contractive  with 
respect  to  the  metric  d  on  X.  If  ^  =  (X; /,^^, .., /^,)  for  some  choice  of  indices 
1  <  Wi  <  W2  <  ...  <  w;  <  iV,  then  we  say  that  ^  is  a  subIFS  of  T . 

Given  an  IFS  Q  =  (X;  gi,  ...,gM)  and  a  sequence  of  indices  p  =  PiP2---Pmi  where 
each  p^  belongs  to  {1,2,  we  can  construct  a  mapping  ^  H  by 

defining 

(BpJ  ,  for  all  B  =  (Bi,B2,  ...,Bv)  e  H^. 

In  a  similar  manner,  given  a  set  of  subIFSs  {^i,...,^l}  of  F,  each  consisting  of 
M  functions,  we  can  construct  mappings  from  to  itself.  Let  a  =  f7i(T2---0'v  € 
{1,  2, ...,  L}'^,  let  p  be  a  E  X  M  matrix  whose  entries  belong  to  [1,2,  and 

here  let  p^  denote  the  row  of  p.  Then  we  define  a  mapping  ^ 
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by 

We  denote  the  sequence  of  all  such  mappings  by  {Hj  :  j  G  J}  where  J  is  the  set  of 
all  indices  (p,cr),  in  some  order.  We  call  =  (H^;  {Hj  :  j  G  J})  the  V -variable 
superlFS  associated  with  the  set  of  subIFSs  of  T . 

We  write  to  denote  the  component  of  i?  G  If  the  space  is  equipped 
with  the  metric  D{B,  C)  :=  max„  {h{By,  Gy)}  where  h  is  the  Hausdorff  metric  on 
H,  then  is  a  complete  metric  space.  The  following  theorem  summarizes 

basic  information  about  More  information  is  presented  in  [4]  and  [7]. 

Theorem  5.  [7]  Let  ©^^^  denote  the  V— variable  superlFS  (H^;  {Hj  :  j  G  J}). 

(i)  If  the  underlying  IFS  T  is  contractive  then  the  IFS  contractive. 

(a)  The  unique  attractor  I  G  H  of  >3^^^  consists  of  a  set  of  V -tuples 

of  compact  subsets  of  X,  and  :=  {By  :  B  G  A^'^\v  =  1,2,...,M}  =  {By  : 

B  G  for  all  v  =  1,2,  ...,y.  (Symmetry  of  the  superlFS  with  respect  to  the 

V  coordinates  causes  this.)  Each  element  of  is  a  union  of  transformations, 
belonging  to  T ,  of  at  most  V  other  elements  of  A^^K 

(Hi)  7/{Afc}^Q  denotes  a  random  orbit  of  Aq  G  under  and  G  El  de¬ 
notes  the  first  component  of  Ak,  then  (with  probability  one)  lim/^^oo  {Ak\  = 

A^^'>  where  the  limit  is  taken  with  respect  to  the  Hausdorff  metric  on  EI(EI(X)). 

Statement  (i)  implies  that  ©^^^  possesses  a  unique  attractor  A^^^  and  that  we 
can  describe  it  in  terms  of  the  chaos  game.  This  technique  is  straightfoward  to 
apply,  since  we  need  only  to  select  independently  at  each  step,  the  indices  p  and  a; 
the  functions  themselves  are  readily  built  up  from  those  of  the  underlying  IFS  T . 

Statement  (ii)  tells  us  that  it  is  useful  to  focus  on  the  set  A^^^  of  first  components 
of  elements  of  It  also  implies  that,  given  any  A  G  and  any  positive  integer 

K,  there  exist  Ai,A2,...,Ay  G  A^^^  such  that  A  =  U;g{y2,...,v}  Ucrco  /o-|ic(A)- 
That  is,  at  any  depth  AT,  A  is  a  union  of  contractions  applied  to  V  sets,  all  belonging 
to  at  most  V  of  which  are  distinct.  In  view  this  property,  the  elements  of 

A^^^  are  called  V-variable  fractal  sets  and  we  refer  to  A^^^  itself  as  a  superfractal. 

Statement  (in)  tells  us  that  we  can  use  random  orbits  {A^j^p  of  Aq  G  EI^  under 
to  sample  the  superfractal  A^^K 

The  idea  of  address  structures,  tops  and  fractal  transformations  can  be  extended 
to  the  individual  sets  that  comprise  A^'^^  see  [4].  We  are  thus  able  to  render  colorful 
images  of  sequences  of  elements  of  A^^^  generated  by  a  more  elaborate  chaos  game 
involving,  at  each  step  a  (M  +  l)-tuple  of  sets,  one  of  which  is  used  to  define  the 
picture  whose  colors  are  used  to  render  the  other  sets. 

Figure  12  illustrates  elements  taken  from  such  an  orbit.  In  this  case  a  2-variable 
superlFS  is  used:  it  consists  of  two  subIFS  of  a  projective  IFS  F,  consisting  of  five 
functions,  detailed  in  [4] .  The  fern-like  structure  of  all  elements  of  the  corresponding 
superfractal  is  ensured  by  a  generalized  version  of  the  collage  theorem. 

7.2.  The  problem  solved  by  superfractals.  The  diverse  forms  that  illustrate  a 
successful  idea,  the  plants  of  a  creative  system,  may  change  as  time  goes  forward. 
The  ideas  of  an  earlier  era  of  geometry,  popular  in  applications  included  cissoids, 
strophoids,  nephroids,  and  astroids:  more  recently  you  would  hear  about  manifolds, 
Ricci  curvature,  and  vector  bundles;  today  you  are  just  as  likely  to  hear  about 
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Figure  12. 
Figure  5. 


Elements  of  a  2-variable  superfractal,  rendered  as  in 


fractals.  Why?  As  technology  advances,  some  applications  become  extinct,  and 
new  ones  emerge. 

In  The  Fractal  Geometry  of  Nature,  [24] ,  Mandelbrot  argues  that  random  frac¬ 
tals  provide  geometrical  models  for  naturally  occurring  shapes  and  forms,  such  as 
coastlines,  clouds,  lungs,  trees,  and  Brownian  motion.  A  random  fractal  is  a  sta¬ 
tistically  self-similar  object  with  non-integer  Hausdorff  dimension.  Although  there 
are  mathematical  theories  for  families  of  random  fractals,  see  for  example  [25],  they 
are  generally  cumbersome  to  use  in  geometric  modelling  applications. 

For  example,  consider  the  problem  of  modelling  real  ferns:  ferns  look  different  at 
different  levels  of  magnification,  and  the  locations  of  the  fronds  are  not  according 
to  some  strictly  deterministic  pattern,  as  in  a  geometric  series,  but  rather  they  have 
elements  of  randomness.  It  seems  that  a  top-down  hierarchical  description,  starting 
at  the  coarsest  level,  and  working  down  to  finer  scales,  is  needed  to  provide  specific 
geometrical  information  about  structure  at  all  levels  of  magnification.  This  presents 
a  problem:  clearly  it  is  time  consuming  and  expensive  in  terms  of  the  amount  of 
data  needed,  to  describe  even  a  single  sample  from  some  statistical  ensemble  of 
such  objects. 

Superfractals  solve  this  problem  by  restricting  the  type  of  randomness  to  be  In¬ 
variable.  This  approach  enables  the  generalized  chaos  game  algorithm,  described 
above,  to  work,  yielding  sequences  of  samples  from  a  probability  distribution  on 
E- variable  sets  belonging  to  a  superfractal.  In  turn,  this  means  that  we  can  approx¬ 
imate  fully  random  fractals  because,  in  the  limit  as  V  tend  to  infinity,  E-variable 
fractals  become  random  fractals  in  the  sense  of  [25],  see  for  example  [7],  Theorem 
51. 

Thus,  we  are  able  to  compute  arbitrarily  accurate  sequences  of  samples  of  ran¬ 
dom  fractals.  Furthermore  we  have  modelling  tools,  obtained  by  generalizing  those 
that  belong  to  IFS  theory,  such  as  collage  theorem  and  fractal  transformations, 
which  extend  in  natural  ways  to  the  E-variable  setting.  In  some  cases  the  Haus¬ 
dorff  dimension  of  these  objects  can  also  be  specified  as  part  of  the  model.  This 
provides  an  approach  to  modelling  many  naturally  occurring  structures  that  is  both 
mathematically  satisfying  and  computationally  workable.  In  particular,  we  see  how 
the  IFS  creative  system  begat  a  new,  even  more  powerful  system,  with  diverse  po¬ 
tential  applications.  This  completes  my  argument  that  iterated  function  systems 
comprise  a  creative  system. 
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8.  Further  reading 

I  would  have  liked  to  tell  you  much  more  about  IFS  theory.  But  this  is  not  a 
review  article,  even  of  some  of  my  own  work.  It  does  not  touch  the  full  range  of  the 
subject,  let  alone  the  mathematics  of  fractal  geometry  as  a  whole.  The  contents 
were  chosen  primarily  to  illustrate  the  idea  of  a  mathematical  creative  system. 

To  survey  mathematical  fractal  geometry,  I  mention  the  series  of  four  conference 
proceedings  [36],  [37],  [38],  and  [39],  carefully  edited  by  Christoph  Bandt,  Martina 
Zahle  and  others.  The  books  by  Falconer,  for  example  [14]  and  [15],  are  good 
textbooks  for  core  material.  A  recent  development  has  been  the  discovery  of  how 
to  construct  harmonic  functions  and  a  calculus  on  certain  fractal  sets,  see  [20]. 
This  was  reported  in  the  Notices  [33].  A  light  introduction  is  [34].  Fractals  and 
number  theory  is  an  important  area;  see  for  example  [22],  [11],  and  [23].  The  topic 
of  noncommutative  fractal  geometry  is  another  fascinating  new  area  [21].  Fractal 
geometry  is  rich  with  creative  possibilities. 

9.  Conclusion 

In  this  essay  I  have  illustrated  the  notion  that  mathematical  ideas  that  survive 
are  creative  systems  in  their  own  right,  with  attributes  that  parallel  some  of  natural 
evolution. 

Creative  systems  dehne,  via  their  DNA,  diverse  forms  and  structures.  There  are 
three  concepts  here:  seeds,  plants,  and  diversity.  Individual  plants  are  products  of 
the  system,  representatives  of  its  current  state  and  utility.  The  system  itself  may 
remain  constant  but  the  plants  evolve,  adapting  to  new  generations  of  minds.  The 
IFS  creative  system  lives  in  my  mind.  But  mainly  I  watch  its  plants:  ones  that 
preoccupy  me  now  are  not  the  same  as  the  ones  that  I  looked  at  years  ago;  the 
crucial  element  is  the  creative  system,  not  the  fractal  fern. 

Plants  provide  the  first  wave  of  conquest  of  new  environments;  an  adapted  version 
of  the  underlying  new  idea  may  follow  later.  The  diversity  of  plants  suggests  a 
multitude  of  possibilities.  Their  seeds  get  into  the  minds  of  engineers  and  scientists. 
Later,  the  underlying  idea,  the  creative  system  itself,  may  take  hold. 

I  think  of  a  good  mathematical  mind,  a  strong  mathematics  department,  and 
a  successful  conference,  as  each  being  like  Darwin’s  bank,  a  rich  ecosystem,  a  fer¬ 
tile  environment  where  ideas  interact  and  diverse  species  of  plants  are  in  evidence. 
Some  of  these  plants  may  be  highly  visible  because  they  can  be  represented  using 
computer  graphics,  while  others  are  more  hidden:  you  may  only  see  them  in  collo- 
quia,  a  few  glittering  words  that  capture  and  describe  something  wonderful,  jump 
from  brain  to  brain  and  there  take  root.  (I  think  of  the  first  time  I  heard  about 
the  Propp-Wilson  algorithm.) 

A  good  mathematical  idea  is  a  creative  system,  a  source  of  new  ideas,  as  rich  in 
their  own  right  as  the  original.  The  idea  that  survives  is  one  that  takes  root  in  the 
minds  of  others:  it  does  so  because  it  is  accessible  and  empowering.  Such  an  idea 
is  likely  to  lead  to  applications,  but  this  applicability  is  more  a  symptom  that  the 
idea  is  a  creative  system,  rather  than  being  causative.  A  good  idea  allows,  invites, 
surprises,  simplihes,  and  shares  itself  without  ever  becoming  smaller;  as  generous, 
mysterious,  and  bountiful  as  nature  itself. 
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Abstract.  The  two  main  theorems  of  this  paper  provide  a 
characterization  of  hyperbolic  affine  iterated  function  systems 
defined  on  R™'.  Atsushi  Kameyama  {Distances  on  Topological 
Self-Similar  Sets,  Proceedings  of  Symposia  in  Pure  Mathe¬ 
matics,  Volume  72.1,  2004)  asked  the  following  fundamental 
question:  given  a  topological  self-similar  set,  does  there  exist 
an  associated  system  of  contraction  mappings?  Our  theorems 
imply  an  affirmative  answer  to  Kameyama’s  question  for  self¬ 
similar  sets  derived  from  affine  transformations  on  R"*. 


1.  Introduction 

The  goal  of  this  paper  is  to  prove  and  explain  two  theorems 
that  classify  hyperbolic  affine  iterated  function  systems  defined  on 
M™.  One  motivation  was  the  question:  when  are  the  functions  of 
an  affine  iterated  function  systems  (IFS)  on  contractions  with 
respect  to  a  metric  equivalent  to  the  usual  euclidean  metric? 

Theorem  1.1  (Classification  for  Affine  Hyperbolic  IFSs).  If  J-  = 
(M™'; /i, /2, /at)  is  an  affine  iterated  function  system,  then  the 
following  statements  are  equivalent. 

(1)  T  is  hyperbolic. 

(2)  T  is  point-fibred. 

(3)  fF  has  an  attractor. 

(4)  T  is  a  topological  contraction  with  respect  to  some  convex 
body  A  C  M™. 
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(5)  T  is  non-antipodal  with  respect  to  some  convex  body  K  C 


Statement  (1)  is  a  metric  condition  on  an  affine  IFS,  statements 
(2)  and  (3)  are  in  terms  of  convergence,  and  statements  (4)  and  (5) 
are  in  terms  of  concepts  from  convex  geometry.  The  terms  con¬ 
tractive,  hyperbolic,  point-hbred,  attractor,  topological  contrac¬ 
tion,  and  non-antipodal  are  defined  in  Definitions  12.21  12.31  12.51 
12.71  15.81  16.51  respectively.  This  theorem  draws  together  some  of 
the  main  concepts  in  the  theory  of  iterated  function  systems.  Ba¬ 
nach’s  classical  Contraction  Mapping  Theorem  states  that  a  con¬ 
traction  /  on  a  complete  metric  space  has  a  fixed  point  xq  and  that 
xq  =  limk^ao  f°^(x),  independent  of  x,  where  ok  denotes  the 
iteration.  The  notion  of  hyperbolic  generalizes  to  the  case  of  an  IFS 
the  contraction  property,  namely  an  IFS  is  hyperbolic  if  there  is  a 
metric  on  M™,  equivalent  to  the  usual  one,  such  that  each  function 
in  the  IFS  is  a  contraction.  The  notion  of  point-hbred,  introduced 
by  Kieninger  [9] ,  is  the  natural  generalization  of  the  limit  condition 
above  to  the  case  of  an  IFS.  While  traditional  discussions  of  fractal 
geometry  focus  on  the  existence  of  an  attractor  for  a  hyperbolic 
IFS,  Theorem  11.11  establishes  that  the  more  geometrical  (and  non¬ 
metric)  assumptions  -  topologically  contractive  and  non-antipodal 
-  can  also  be  used  to  guarantee  the  existence  of  an  attractor.  Ba¬ 
sically  a  function  /  :  M™'  — >  is  non-antipodal  if  certain  pairs 

of  points  (antipodal  points)  on  the  boundary  of  K  are  not  mapped 
by  /  to  another  pair  of  antipodal  points. 

Since  the  implication  (1)  ^  (2)  is  the  Contraction  Mapping  The¬ 
orem  when  the  IFS  contains  only  one  affine  mapping,  Theoremll.il 
contains  an  affine  IFS  version  of  the  converse  to  the  Contraction 
Mapping  Theorem.  Thus,  our  theorem  provides  a  generalization  of 
results  proved  by  L.  Janos  [7]  and  S.  Leader  m-  Such  a  converse 
statement  in  the  IFS  setting  has  remained  unclear  until  now. 

Although  not  every  affine  IFS  JF  =  (M"^;/i,/2,  .•■,/iv)  is  hyper¬ 
bolic  on  all  of  M™,  the  second  main  result  states  that  if  T  has  a 
coding  map  (Definition  12. 4p.  then  T  is  always  hyperbolic  on  some 
affine  subspace  of  M™. 

Theorem  1.2.  If  T  =  (M™'; /i, /2, ..., /w)  o,n  affine  IFS  with  a 
coding  map  tt  :  D  ^  M"*,  then  J-  is  hyperbolic  on  the  affine  hull  of 
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7r(r2).  In  particular,  if  contains  a  non-empty  open  subset  of 

M™,  then  T  is  hyperbolic  on  M™'. 

Although  he  used  slightly  different  terminology,  Kameyama  [8] 
posed  the  following  FUNDAMENTAL  QUESTION:  Is  an  affine 
IFS  with  a  coding  map  vr  :  — >  MN  hyperbolic  when  restricted  to 
7r(0)  ?  An  affirmative  answer  to  this  question  follows  immediately 
from  Theorem  11.21 

Our  original  motivation,  however,  was  not  Kameyama’s  question, 
but  rather  a  desire  to  approximate  a  compact  subset  T  C  M™'  as  the 
attractor  A  of  an  iterated  function  system  T  =  (M™; /i, /2,  .../at), 
where  each  /„  :  — >  MN  is  affine.  This  task  is  usually  done  using 

the  “collage  theorem”  m  by  choosing  an  IFS  T  so  that  the 
Hausdorff  distance  dm{T,  T  (T))  is  small.  If  the  IFS  T  is  hyperbolic, 
then  we  can  guarantee  it  has  an  attractor  A  such  that  d^{T,A)  is 
comparably  small.  But  then  the  question  arises:  how  does  one 
know  if  T  is  hyperbolic? 

The  paper  is  organized  as  follows.  Section  2  contains  notation, 
terminology,  and  definitions  that  will  be  used  throughout  the  pa¬ 
per.  Section  3  contains  examples  and  remarks  relating  iterated 
function  systems  and  their  attractors  to  Theorems  11.11  and  11.21  In 
Example  3.1  we  show  that  an  affine  IFS  can  be  point-fibred,  but 
not  contractive  under  the  usual  metric  on  M™.  Thus,  some  kind 
of  remetrization  is  required  for  the  system  to  be  contractive.  In 
Example  3.2  we  show  that  an  affine  IFS  can  contain  two  linear 
maps  each  with  real  eigenvalues  all  with  magnitudes  less  than  1, 
but  still  may  not  be  point-hbered.  Thus,  Theorem  11.11  cannot  be 
phrased  only  in  terms  of  eigenvalues  and  eigenvectors  of  the  indi¬ 
vidual  functions  in  the  IFS.  Indeed,  in  Example  3.3  we  explain  how, 
given  any  integer  M  >  0,  there  exists  a  linear  lES  (M^;  Li,  L2)  such 
that  each  operator  of  the  form  with  cjj  G  {1)2}  for 

j  =  l,2,...,/c,  and  k  <  M,  has  spectal  radius  less  than  one,  while 
LiL^  has  spectral  radius  larger  than  one.  This  is  related  to  the 
joint  spectral  radius  [16]  of  the  pair  of  linear  operators  Li,L2  and 
to  the  associated  hniteness  conjecture,  see  for  example  m-  In  Sec¬ 
tion  8  we  comment  on  the  relationship  between  the  present  work 
and  recent  results  concerning  the  joint  spectral  radius  of  finite  sets 
of  linear  operators.  Example  3.4  provides  an  affine  IFS  on  that 
has  a  coding  map  tt,  but  is  not  point-fibred  on  and  hence  by 
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Theoremll.il  not  hyperbolic  on  It  is,  however,  point-fibred  and 
hyperbolic  when  restricted  to  the  x-axis,  which  is  the  affine  hull  of 
7r(II),  thus  illustrating  Theorem II .21 

For  the  proof  of  Theorem  1 1.1 1  we  provide  the  following  roadmap. 

(1)  The  proof  that  statement  (1)  statement  (2)  is  provided 
in  Theorem  14.11 

(2)  The  proof  that  statement  (2)  statement  (3)  is  provided 
in  Theorem  14.31 

(3)  The  proof  that  statement  (3)  statement  (4)  is  provided 
in  Theorem  15.101 

(4)  The  proof  that  statement  (4)  statement  (5)  is  provided 
in  Proposition  16.61 

(5)  The  proof  that  statement  (5)  statement  (1)  is  provided 
in  Theorem  16.71 

Theorem  0  is  proved  in  section  7. 


2.  Notation  and  Definitions 

We  treat  M”*  as  a  vector  space,  an  affine  space,  and  a  metric 
space.  We  identify  a  point  x  =  {xi,X2,  with  the 

vector  whose  coordinates  are  xi,X2,  ■■■,Xm-  We  write  0  G  for 
the  point  in  whose  coordinates  are  all  zero.  The  standard  basis 
is  denoted  {ei ,  62,  ■  ■  ■ ,  6^}-  The  inner  product  between  x,  y  G  R™  is 
denoted  by  (x,  y).  The  2-norm  of  a  point  x  G  R"*  is  ||x||2  =  (x,  x), 

and  the  euclidean  metric  dg  :  R™'  x  R™  — >  [0,oo)  is  defined  by 
d_B(x,y)  =  ||x  — y||2  forallx,y  G  R”*.  The  following  notations, 
conventions,  and  definitions  will  also  be  used  throughout  this  paper: 

(1)  A  convex  body  is  a  compact  convex  subset  of  R™  with  non¬ 
empty  interior. 

(2)  For  a  set  B  in  R™,  the  notation  conv{B)  is  used  to  denote 
the  convex  hull  of  B. 

(3)  For  a  set  B  G  R™',  the  affine  hull,  denoted  aff(i?),  of  B  is  the 
smallest  affine  subspace  containing  B,  i.e.,  the  intersection 
of  all  affine  subspaces  containing  B. 

(4)  The  symbol  BI  will  denote  the  nonempty  compact  subsets 
of  R™,  and  the  symbol  de  will  denote  the  Hausdorff  metric 
on  EI.  Recall  that  (R”*,de)  is  a  complete  metric  space. 
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(5)  A  metric  d  on  M™'  is  said  to  be  Lipschitz  equivalent  to  dE 
if  there  are  positive  constants  r  and  R  such  that 

rdE{x,y)  <  d{x,y)  <  RdE{x,y), 

for  all  x,y  €  M™'.  If  two  metrics  are  Lipschitz  equivalent, 
then  they  induce  the  same  topology  on  M™',  but  the  converse 
is  not  necessarily  true. 

(6)  For  any  two  subsets  A  and  B  of  the  notation  A  —  B  := 

{x  —  y  :  x  G  A  and  y  G  5}  is  used  to  denote  the  pointwise 
subtraction  of  elements  in  the  two  sets. 

(7)  For  a  positive  integer  N,  the  symbol  =  {1,  2, . . .  ,  N}°°  will 
denote  the  set  of  all  infinite  sequences  of  symbols  {crk}'^=i 
belonging  to  the  alphabet  {1,2,..., A"}.  The  set  Q,  is  en¬ 
dowed  with  the  product  topology.  An  element  of  a  G  O 
will  also  be  denoted  by  the  concatenation  a  =  (T1CJ2CT3  . . . , 
where  ak  denotes  the  component  of  a.  Recall  that  since 
n  is  endowed  with  the  product  topology,  it  is  a  compact 
Hausdorff  space. 

Definition  2.1  (IFS).  If  A  >  0  is  an  integer  and  fn  ■  K™'  ^  K™', 
n  =  1,  2, . . . ,  A,  are  continuous  mappings,  then  T  =  (M™;  /i,  /2, ...,  /at) 
is  called  an  iterated  function  system  (IFS).  If  each  of  the  functions 
in  T  is  an  affine  map  on  then  T  is  called  an  affine  IFS. 

Definition  2.2  (Contractive  IFS).  An  IFS  T  =  (M™';  /i,  /2, ...,  /n) 
is  contractive  when  each  /„  is  a  contraction.  Namely,  there  is  a 
number  G  [0,1)  such  that  dE{fn{x),  fn{y))  <  CindE{x,y)  for  all 
x,y  €  M™',  for  all  n. 

Definition  2.3  (Hyperbolic  IFS).  An  IFS  T  =  (M™';/i,/2,  ...,/Ar) 
is  called  hyperbolic  if  there  is  a  metric  on  Lipschitz  equivalent 
to  the  given  metric  so  that  each  /„  is  a  contraction. 

Definition  2.4  (Coding  Map).  A  continuous  map  n  :  H  ^  is 
called  a  coding  map  for  the  IFS  T  =  (R™^;  /i,  /2,  •••,  /at)  if,  for  each 
n  =  1,  2, . . . ,  A,  the  following  diagram  commutes, 

(2.1)  TT  i  i  TT 

fn 

where  the  symbol  :  II  ^  11  denotes  the  inverse  shift  map  defined 
by  Sn{cr)  =  na. 
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The  notion  of  a  coding  map  is  due  to  J.  Kigami  [10]  and  A. 
Kameyama  [8]. 

Definition  2.5  (Point-Fibred  IFS).  AnIFS  T  =  (M™';  /i,  /2, /at) 
is  point-fibred  if,  for  each  cr  =  cti  cJ2  CTs  •  •  •  G  the  limit  on  the  right 
hand  side  of 

(2.2)  7r{a)  :=  lim 

k^OQ 

exists,  is  independent  of  x  G  M™'  for  fixed  a,  and  the  map  vr  :  — > 

is  a  coding  map. 

It  is  not  difficult  to  show  that  12.21  is  the  unique  coding  map  of 
a  point-fibred  IFS.  Our  notion  of  a  point-fibred  iterated  function 
system  is  similar  to  Kieninger’s  Definition  4.3.6  [9],  p.97.  However, 
we  work  in  the  setting  of  complete  metric  spaces  whereas  Kieninger 
frames  his  definition  in  a  compact  Hausdorff  space. 

Definition  2.6  (The  Symbol  J^{B)  for  an  IFS).  For  an  IFS  T  = 
(M'";  /i,  /2, ...,  /iv)  define 

N 

HB)  =  U  fn{B). 

n=\ 

(We  use  the  same  symbol  J-  both  for  the  IFS  and  the  mapping.) 
For  H  G  H,  let  denote  the  fc-fold  composition  of  J-,  i.e.,  the 

union  of  o  o  ■  ■  ■  o  over  all  words  cricr2  •  •  •  cr^  of  length 

k. 

Definition  2.7  (Attractor  for  an  IFS).  A  set  A  G  El  is  called  an 
attractor  of  an  IFS  T  =  (M™;  fi,  f2,  ■■■,  In)  if 

(2.3)  A  =  B(A) 
and 

(2.4)  A  =  lim 

k^oo 

the  limit  with  respect  to  the  Hausdorff  metric,  for  all  H  G  El. 

If  an  IFS  has  an  attractor  A,  then  clearly  A  is  the  unique  attrac¬ 
tor.  It  is  well  known  that  a  hyperbolic  IFS  has  an  attractor.  An 
elegant  proof  of  this  fact  is  given  by  J.  Hutchinson  [6].  He  observes 
that  a  contractive  IFS  B  induces  a  contraction  T  :  El  — >  El,  from 
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which  the  result  follows  by  the  contraction  mapping  theorem.  See 
also  M.  Hata  [5]  and  R.  F.  Williams  m- 

In  section  4  it  is  shown  that  a  point-fibred  IFS  T  has  an  attractor 
A,  and,  moreover,  if  vr  is  the  coding  map  of  J- ^  then  A  =  7r(n). 
Often  o  is  considered  as  the  “address”  of  the  point  vr(cj)  in  the 
attractor.  In  the  literature  on  fractals  (for  example  J.  Kigami  [TO]) 
there  is  an  approach  to  the  concept  of  a  self-similar  system  without 
reference  to  the  ambient  space.  This  approach  begins  with  the  idea 
of  a  continuous  coding  map  tt  and,  in  effect,  defines  the  attractor 
as  tt{Q). 

3.  Examples  and  Remarks  on  Iterated  Function  Systems 

This  section  contains  examples  and  remarks  relevant  to  Theo¬ 
rems  11.11  and  O 


EXAMPLE  3.1  [A  Point-fibred,  not  Contractive  IFS]  Consider 
the  affine  IFS  consisting  of  a  single  linear  function  on  given  by 
the  matrix 


/ 


Note  that  the  eigenvalues  of  /  equal  ±i.  Since 


lim  = 


lim  T 

n^oo 


-1 


ar 


2 

0 

i-hr 


T  = 


0  0 
0  0 


where  T  is  the  change  of  basis  matrix,  this  IFS  is  point-fibred. 
However,  since 


the  mapping  is  not  a  contraction  under  the  usual  metric  on  M^. 
Theorem  m  however,  guarantees  we  can  remetrize  R^  with  an 
equivalent  metric  so  that  /  is  a  contraction. 


EXAMPLE  3.2  [An  IFS  with  Point-Fibred  Functions  that  is  not 
Point-Fibred]  In  the  literature  on  affine  iterated  function  systems, 
it  is  sometimes  assumed  that  the  eigenvalues  of  the  linear  parts 
of  the  affine  functions  are  less  than  1  in  modulus.  Unfortunately, 
this  assumption  is  not  sufficient  to  imply  any  of  the  five  statements 
given  in  Theorem  11.11  While  the  affine  IFS  (R'”;/)  is  point-fibred 
if  and  only  if  the  eigenvalues  of  the  linear  part  of  /  all  have  moduli 
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strictly  less  than  1,  an  analogous  statement  cannot  be  made  if  the 
number  of  functions  in  the  IFS  is  larger  than  1. 

Consider  the  affine  IFS  T  =  (M^;  /i,  72)  ,  where 

/.  =  (?  and  /,=  (“  5). 

As  noted  in  Example  3.1 

lim  /i°”u  =  lim  f 

n— >00  n— >00  \  U  j 

for  any  vector  u.  Thus,  both  T\  =  (M^;/i)  and  T2  =  (M^;/2)  are 
point-fibred.  Unfortunately,  their  product  is  the  matrix 

f\°h  =  ,  so  that  J.^(/io/2)°"  ^  =  -hoo. 

Thus,  the  IFS  T  =  (M^;  /i,  72)  fails  to  be  point-fibred. 

Remark  3.1. 

While  it  is  true  that  (1)  (2)  in  Theorem  11.11  even  without 

the  assumption  that  the  IFS  is  affine,  the  converse  is  not  true  in 
general.  Kameyama  [8]  has  shown  that  there  exists  a  point-fibred 
IFS  that  is  not  hyperbolic.  We  next  give  an  example  of  an  affine 
IFS  with  a  coding  map  that  is  not  point-fibred.  Thus,  the  set  of 
IFSs  (with  a  coding  map)  strictly  contains  the  set  of  point-fibred 
IFSs  which,  in  turn,  strictly  contains  the  set  of  hyperbolic  IFSs. 


EXAMPLE  3.3  [The  Eailure  of  a  Einite  Eigenvalue  Test  to  Imply 
Point-Eibred]  Consider  the  linear  lES  T  =  (M^;Li,L2)  ,  where 


u  = 


0  2 


and 


r  —  ( ^  ^  _  D 

Usin^  acos0  )  “ 


where  Rq  denotes  rotation  by  angle  9,  and  0  <  a  <  1.  Then  L”  has 
eigenvalues  ±1/2"'  while  the  eigenvalues  of  ±2  both  have  magnitude 
a"  <  1.  Eor  example,  if  we  choose  9  =  vr/S  and  a  =  31/32  then  it  is 
readily  verified  that  the  eigenvalues  of  L1L2  and  L2L1  are  smaller 
than  one  in  magnitude  and  that  one  of  the  eigenvalues  of  L1L2L2 
is  1.4014...  .  Hence,  in  this  case,  the  magnitudes  of  the  eigenvalues 
of  the  linear  operators  Li,  L2,  Tf,  L1L2,  L2L1,  are  all  less  than 
one,  but  ||(LiL2A2)’^  a:||  does  not  converge  when  x  G  is  any 
eigenvector  of  L1L2L2  corresponding  to  the  eigenvalue  1.4014....  It 
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follows  that  the  IFS  Li,  L2)  is  not  point-hbred.  By  using  the 
same  underlying  idea  it  is  straightforward  to  prove  that,  given  any 
positive  integer  M,  we  can  choose  a  close  to  1  and  9  close  to  0  in 
such  a  way  that  the  eigenvalues  of  (where  dj  G  {1,  2} 

for  j  =  1,  2, k,  with  k  <  M)  are  all  of  magnitude  less  than  one, 
while  LiL^  has  an  eigenvalue  of  magnitude  larger  than  one. 

EXAMPLE  3.4  [A  non-Hyperbolic  Affine  IFS]  Let  T  =  (M^;  /g,  /^)  , 
where 

/o(xi,X2)  =  (^Xi,X2),  /i(xi,  X2)  =  {^xi  +  X2). 

This  IFS  has  a  coding  map  tt  with  =  {0, 1}°°  and  vr(cr)  =  (O.cr,  0), 
where  O.cr  is  considered  as  a  base  2  decimal.  Since  limfc^oo  /^-j  o 
/0-2  o  ■  ■  ■  o  /o-j,  (xi,  X2)  =  (O.iT,  X2)  depends  on  the  choice  of  the  points 
(xi,X2)  G  this  IFS  cannot  be  point-fibred.  Hence,  by  Theo¬ 
rem  0  the  IFS  J-  is  also  not  hyperbolic.  However,  it  is  clearly 
hyperbolic  when  restricted  to  the  x-axis,  the  affine  hull  of  unit  inter¬ 
val  7r(H)  =  [0, 1]  X  {0}.  Thus,  this  example  illustrates  TheoremlI.21 

A  key  fact  used  in  the  proof  of  Theorem II. II is  that  the  set  of  an¬ 
tipodal  points  in  a  convex  body  equals  the  set  of  diametric  points. 
The  definitions  of  antipodal  and  diametric  points  are  given  in  Def¬ 
initions  ED  and  16.21  respectively.  The  equality  between  these  two 
point  sets  is  proved  in  Theorem  16.41  While  it  is  possible  that  this 
result  is  present  in  the  convex  geometry  literature,  it  does  not  seem 
to  be  well-known.  For  example,  it  is  not  mentioned  in  the  works  of 
Moszynska  m  or  Schneider  m-  This  equivalence  between  antipo¬ 
dal  and  diametric  points  is  crucial  to  our  work  because  it  provides 
the  remetrization  technique  at  the  heart  of  Theorem  16.71  which 
implies  that  a  non-antipodal  IFS  is  hyperbolic.  A  consequence  of 
Theorem  11.11  is  that  a  non-antipodal  affine  IFS  has  the  seemingly 
stronger  property  of  being  topologically  contractive. 

4.  Hyperbolic  Implies  Point-fibred  Implies  The 
Existence  of  an  Attractor 

The  implications  (1)  =>  (2)  (3)  in  Theorem  11.11  are  proved 

in  this  section.  For  this  section  we  also  introduce  the  notation 
/f7|fc  =  /f7i  o  /crj  o  •  •  •  o  fa^{x).  Note  that,  for  k  fixed,  /^|fc(x)  is  a 
function  of  both  x  and  cr. 
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Theorem  4.1.  If  =  (M"*;/i,/2,  ...j/a)  is  a  hyperbolic  IFS,  then 
T  is  point  fibred. 

Proof.  For  a  G  12,  the  proof  that  the  limit  limfc^oo  fa\k  exists  and  is 
independent  of  x  is  virtually  identical  to  the  proof  of  the  classical 
Contraction  Mapping  Theorem.  Moreover,  the  same  proof  shows 
that  the  limit  is  uniform  in  a. 

With  TT  :  12  — >  M™'  defined  by  vr((T)  =  lim^^oo  fa\k  h  is  easy  to 
check  that,  for  each  n  =  1,  2, . . .  ,  A^,  the  diagram  12.11  commutes. 

It  only  remains  to  show  that  vr  is  continuous.  With  x  fixed, 
fa\kix)  is  a  continuous  function  of  a.  This  is  simply  because,  if 
cr,  r  G  12  are  sufficiently  close  in  the  product  topology,  then  they 
agree  on  the  first  k  components.  By  Definition  12. the  function  vr 
is  then  the  uniform  limit  of  continuous  (in  a)  functions  defined  on 
the  compact  set  12.  Therefore,  vr  is  continuous.  □ 

Let  tF  he  a.  point-fibred  affine  IFS,  and  let  A  denote  the  set 

A  :=  7r(Q). 

According  to  Theorem  14.31  A  is  the  attractor  of  tF. 

Lemma  4.2.  Let  IF  =  (M™; /i, /2, ..., /a)  be  a  point-fibred  affine 
IFS  with  coding  map  tt  :  12  — >  R™' .  If  B  C  R™'  is  compact,  then  the 
convergence  in  the  limit 

7r(cj)  =  lim  /^|fc(x) 
fc— >oo 

is  uniform  in  u  =  criiT2  •  •  •  G  12  and  x  G  B  simultaneously. 

Proof.  Only  the  uniformity  requiress  proof.  Express  fn{x)  =  L„x  + 
On,  where  Ln  is  the  linear  part.  Then 

(4.1) 

fa\k{x)  —  -Lo-jfc(x)  -|-  +  Lo-|fc_2(ao-j,_j)  L„\ia2  +  oi 

-i-^a\k(.X^  4“  f(j\k(.^^‘ 

From  equation  14.11  it  follows  that,  for  any  x,y  €  B, 

(4.2) 

dE{fa\k{x),  fa\k{y))  =  \\La\kix  -  y)\\.^ 

m 

<  sup  (E  2  \cj\  ||Lo-|fc(ej)||2  :  cici 

J=1 

<  c  max  \\fa\k{ej)  -  /f7|fc(0)||2  , 
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where  c  =  2m  ■  sup {maxj  \cj\  :  ciei  +  •  ■  ■  +  CmCm  ^  B}  and  where 
is  a  basis  for  M’". 

Let  e  >  0.  From  the  definition  of  point-fibred  there  is  a  fcj, 
independent  of  a,  such  that  if  fc  >  kj,  then 

\\fa\k{ej)  -  7r(cT)||2  <  ^  and  \\fa\ki0)  -  7r(o-)||2  < 

which  implies  ||/(T|fc(ej)  — /o-|fc(0)||2  <  This  and  equation  14.21 
implies  that  if  A;  >  A:  :=  maxj  kj,  then  for  any  x,y  &  B  we  have 

(4.3)  dEiL\k{x),f^\kiy))  <  =  |- 

Let  6  be  a  fixed  element  of  B.  There  is  a  kb,  independent  of  a,  such 
that  if  A:  >  kb,  then  dEifa\k{b),iT{a))  <  |.  If  fe  >  max{kb,k)  then, 
by  equation  14. 3t  for  any  x  G  B 

dE{fa\kix),TT{a))  <  dEifa\k{x),fa\k{b))  +  dEifa\k{b),E{a))  <  ^  +  ^  = 

□ 


Theorem  4.3  (A  Point-Fibred  IFS  has  an  Attractor).  If  B  = 
(M™'; /i, /2, /at)  is  a  point-fibred  affine  IFS,  then  T  has  an  at- 
traetor  A  =  vr(n),  where  tt  :  If  — >  M™’  is  the  eoding  map  of  IF . 

Proof.  It  follows  directly  from  the  commutative  diagram  (I2.1h  that 
A  obeys  the  self-referential  equation  (12  ■3h.  We  next  show  that  A 
satisfies  equation  (12. 4h. 

Let  e  >  0.  We  must  show  that  there  is  an  M  such  that  if 
k  >  M,  then  dM{F°^{B),Tr{Il))  <  e.  It  is  sufficient  to  let  M  = 
max(Mi,M2),  where  Mi  and  M2  are  defined  as  follows. 

First,  let  a  be  an  arbitrary  element  of  A.  Then  there  exists  a 
cr  G  n  such  that  a  =  vr(cr).  By  Lemma  14.21  there  is  an  Mi  such 
that  if  A:  >  Mi,  then  dE{fa\k{b),a)  =  dE{fa\k{b),T:{a))  <  e,  for  all 
b  ^  B.  In  other  words,  A  lies  in  an  e-neighborhood  of  F°^{B). 

Second,  let  b  be  an  arbitrary  element  of  B  and  a  an  arbitrary 
element  of  II.  If  a  :=  vr(cj)  G  A,  then  there  is  an  M2  such  that 
if  Ac  >  M2,  then  dE{fa\k{b)F)  =  dE{fa\k{b),'T^ffi))  <  e.  In  other 
words,  F°^{B)  lies  in  an  e- neighborhood  of  A.  □ 
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5.  An  IFS  WITH  AN  Attractor  is  Topologically 
Contractive 

The  goal  of  this  section  is  to  establish  the  implication  (3)  ^  (4) 
in  Theoremll.il  We  will  show  that  if  an  affine  IFS  has  an  attractor 
as  defined  in  Defintion  12.71  ,  then  it  is  a  topological  contraction. 
The  proof  uses  notions  involving  convex  bodies. 

Definition  5.1.  A  convex  body  K  is  centrally  symmetric  if  it  has 
the  property  that  whenever  x  G  K,  then  —x  €  K. 

A  well-known  general  technique  for  creating  centrally  symmetric 
convex  bodies  from  a  given  convex  body  is  provided  by  the  next 
proposition. 

Proposition  5.2.  If  a  set  K  is  a  convex  body  in  M™',  then  the  set 
K'  =  K  —  K  is  a  centrally  symmetric  convex  body  in  M"*. 

Definition  5.3  (Minkowski  Norm).  If  A  is  a  centrally  symmetric 
convex  body  in  M™,  then  the  Minkowski  norm  on  M™'  is  defined  by 

||x||^  =  inf  {A  >  0  :  X  G  A  A}. 

The  next  proposition  is  also  well-known. 

Proposition  5.4.  If  K  is  a  centrally  symmetric  convex  body  in 
M™',  then  the  function  \\x\\j^  defines  a  norm  on  M™'.  Moreover,  the 
set  A  is  the  unit  ball  with  respect  to  the  Minkowski  norm  ||x||^. 

Definition  5.5  (Minkowski  Metric).  If  A  is  a  centrally  symmetric 
convex  body  in  M™’  and  ||x||^  is  the  associated  Minkowski  norm, 
then  define  the  Minkowski  metric  on  M™'  by  the  rule 

dK{x,y)  :=  ||x  -  y||^  . 

While  R.  Rockafeller  [13]  refers  to  such  a  metric  as  a  Minkowski 
metric,  the  reader  should  be  aware  that  this  term  is  also  associated 
with  the  metric  on  space-time  in  theory  of  relativity.  Since,  for  any 
convex  body  A  there  are  positive  numbers  r  and  R  such  that  A 
contains  a  ball  of  radius  r  and  is  contained  in  a  ball  of  radius  R, 
the  following  proposition  is  clear. 

Proposition  5.6.  If  d  is  a  Minkowski  metric,  then  d  is  Lipschitz 
equivalent  to  the  standard  metric  dE  on  M™. 
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Proposition  5.7.  A  metric  d  :  M™’  x  — >  [0,  oo)  is  a  Minkowski 

metric  if  and  only  if  it  is  translation  invariant  and  distances  behave 
linearly  along  line  segments.  More  specifically, 

(5.1) 

d{x  +  z,y  +  z)  =  d{x,y)  and  d{x,{l  —  X)x  +  \y)  =  Xd{x,y) 

for  all  x,y,z  €  M™'  and  all  A  G  [0, 1] . 

Proof.  For  a  proof  see  Rockafeller  [14]  pp. 131-132.  □ 

Definition  5.8  (Topologically  Contractive  IFS).  An  IFS  T  = 

{M"*; /i, /2, ..., /at}  is  called  topologically  contractive  if  there  is  a 
convex  body  K  such  that  T  (K)  C  int{K). 

The  proof  of  Theorem  15.101  relies  on  the  following  lemma  which 
is  easily  proved. 

Lemma  5.9.  If  g  :  ^  is  affine  and  S  C  M™',  then 

g{conv{S))  =  conv{g{S)). 

Theorem  5.10  (The  Existence  of  an  Attractor  Implies  a  Topo¬ 
logical  Contraction).  For  an  affine  IFS  !F  =  {M™';  /i,  f^, ...,  fw},  if 
there  exists  an  attractor  A  G  H  of  the  affine  IFS  T  =  {M™';  /i,  f^,  /at}; 

then  IF  is  topologically  contractive. 

Proof.  The  proof  of  this  theorem  unfolds  in  three  steps. 

(1)  There  exists  a  convex  body  Ki  and  a  positive  integer  t  with 
the  property  that  F°^  (Ai)  C  int  (Ki). 

(2)  The  set  Ki  is  used  to  define  a  convex  body  K2  such  that 
An  (A2)  C  int{K2),  where  fn{x)  =  LnX  +  an  and  n  = 
1,2,... ,iV. 

(3)  There  is  a  positive  constant  c  such  that  the  set  K  =  CK2 
has  the  property  F  (K)  C  int  (A). 

Proof  of  Step  (1).  Let  A  denote  the  attractor  of  F .  Let  Ap  =  {x  G 
:  dH({a^}5A)  <  p}  denote  the  dilation  of  A  by  radius  p  >  0. 
Since  we  are  assuming  hmfc_>oo  d^{F°^[Ap),  A)  =  0,  we  can  find  an 
integer  t  so  that  dM{F°^ {Ai) ,  A)  <  1.  Thus, 


(5.2) 


F°\Ai)  C  int{Ai). 
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If  we  let  Ki  :=  conv  (^i),  then 

{Ki)  =  U  U  •  •  •  U  (/ii  o  o  •  •  •  o  4)  (coni;(^i)) 

=  U  U  ■  "  U  ^-^*1  °  /*2  °  ■  ■  ■  °  fki^i))  (by  Lemma[5JD 

c  u  u-u  conv  {int{Ai))  =  conv{int{Ai))  (by  inclusion  (5.2)) 

C  int{conv{Ai))  =  int  {Ki) . 

This  argument  completes  the  proof  of  Step  (1). 

Proof  of  Step  (2).  Consider  the  set 
t-i 

K2  :=  '^{conv{T°^{Ki)  -  conv{J^°’^{Ki))). 
k=0 

The  set  K2  is  a  centrally  symmetric  convex  body  because  it  is  a 
finite  Minkowski  sum  of  centrally  symmetric  convex  bodies.  If  any 
affine  map  fn  in  T  is  written  fnix)  =  LnX  +  an,  where  Ln  ■ 
denotes  the  linear  part,  then 
t-i 

Ln{K2)  =  Ln  {conv{T°^ {Ki)  —  conv{T°^ {Ki)'^  (since  Ln  is  a  linear  map) 

k={) 
t-1 

=  (conv{Ln  (iiTi)^)  —  conv{Ln  )^  (by  Lemma (5^ 

/c=0 

t-1 

=  yy  (conv{fn  (iPi)^)  —  conv{fn  (since  the  o^s  cancel) 

k=0 
t-1 

C  (^conn(.^°(^+b  (Ki))  -  conv{T°^’^+^\Ki))'^ 

k=0 

t-1 

=  (^conv{L'°^  (Ki)  —  conv{J^°^{Ki))j  +  y^  (conv{L'°^{Ki)  —  conv{!F°^{Ki))^ 

k=l 

t-1 

C  [int  [Ki]  —  int  (i^i))  +  y^^[conv[iF°^ (Kj)  —  conv[iF°^ (Ki)))  (by  Step  1) 

k=l 

=  int[K2). 
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The  second  to  last  inclusion  follows  from  the  fact  that 
fn  {Ki))  C  (Ki).  The  last  equality  follows  from  the 

fact  that  if  O  and  C  are  symmetric  convex  bodies  in  then 
int{0)  +  C  =  int  (O  -h  C).  We  have  now  completed  the  proof  of 
Step  (2). 

Proof  of  Step  (3).  It  follows  from  Step  (2)  and  the  compactness  of 
K2  that  there  is  a  constant  a  G  (0, 1)  such  that  dK2{Ln{x) ■,  Ln{y))  < 
adK2{x,  y)  for  all  x,  y  G  and  all  n  =  1,  2, . . .  , 

Let 

r 


where  r  =  max{(ii^2(oi,  0),  0), . . . ,  dK2  («Af)  0)}.  If  x  G  CK2 

and  /(x)  =  Lx  +  a  is  any  function  in  the  ITS  then 

\\fix)\\K2  =  dK2  (/  (x) ,  0)  =  dK2  {Lx  +  a,  0)  <  dK2  {Lx  +  a,  Lx)  +  dK2  {Lx,  0) 
=  di^2  0)  +  dK2  {Lx,  0)  (by  Equation  (|5.1I)1 
<  r  +  adK2  (x,  0)  =  r  +  a  HxH^^^ 

<r  +  ac<{c  —  ac)  +  ac  =  c. 

This  inequality  shows  that  T  {0X2)  C  int{cK2).  □ 

6.  A  Non- Antipodal  Affine  IFS  is  Hyperbolic 

Let  5”^“^  C  denote  the  unit  sphere  in  M™'.  For  a  convex  body 
K  C  and  u  G  5”^“^  there  exists  a  pair,  {Hu,  ,  of  distinct 
supporting  hyperplanes  of  K,  each  orthogonal  to  u  and  with  the 
property  that  they  both  intersect  dK  but  contain  no  points  of  the 
interior  of  K.  Since  by  definition  a  convex  body  has  non-empty 
interior,  this  pair  will  be  unique.  The  pair  is  usually 

referred  to  as  the  two  supporting  hyperplanes  of  K  orthogonal  to 
u.  (See  Moszynska  [13],  p.l4.) 

Definition  6.1  (Antipodal  Pairs).  If  AT  C  is  a  convex  body 
and  u  G  then  define 

Au  :=  Au{K)  =  {{p,  q)  G  {Hu  n  dK)  x  (H_„  n  BK)}  and 
^  :=  A{K)  =  IJ  A. 

We  say  that  {p,  q)  is  an  antipodal  pair  of  points  with  respect  to  K 
if  {p,  q)  G  A. 
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Definition  6.2  (Diametric  Pairs).  If  X  C  is  a  convex  body, 
and  then  define  the  diameter  of  K  in  the  direction  u  to 

be 

D{u)  =  max{||x  —  y\\2  ■  x,  y^K,  x  —  y  =  au,  a  G  M}. 

The  maximum  is  achieved  at  some  pair  of  points  belonging  to  dK 
because  K  x  K  is  convex  and  compact,  and  ||x  —  1/II2  is  continuous 
for  (x,  y)  ^  K  X  K.  Now  define 

Vu  =  {(p,  q)  G  dK  X  dK  :  D{u)  =  ||g  —  PII2}  and 

V=  IJ  V^. 

We  say  that  (p,  q)  G  is  a  diametric  pair  of  points  in  the  direction 
of  u,  and  that  V  is  the  set  of  diametric  pairs  of  points  of  K. 

Definition  6.3  (Strictly  Convex).  A  convex  body  K  is  strictly 
convex  if,  for  every  two  points  x,y  G  K,  the  open  line  segment 
joining  x  and  y  is  contained  in  the  interior  of  K. 

We  write  xy  to  denote  the  closed  line  segment  with  endpoints  at 
X  and  y  so  that  y  —  x  is  the  vector,  in  the  direction  from  x  to  y, 
whose  magnitude  is  the  length  of  xy. 

Theorem  6.4.  If  K  C  M™'  is  a  convex  body,  then  the  set  of  antipo¬ 
dal  pairs  of  points  of  K  is  the  same  as  the  set  of  diametric  pairs  of 
points  of  K,  i.e., 

A  =  V. 


Proof.  First  we  show  that  A  CV.  If  (p,  q)  G  A,  then  p  G  Hu  H  dK 
and  q  G  H-u  C  dK  for  some  u  G  Clearly  any  chord  of  K 

parallel  to  pq  lies  entirely  in  the  region  between  Hu  and  H-u  and 
therefore  cannot  have  length  greater  than  that  of  pq.  So  D{q—p)  = 
||g  —  pII  and  (p,  g)  G  Hg-p  C  V.  Note,  for  use  later  in  the  proof, 
that  if  K  is  strictly  convex,  then  pq  is  the  unique  chord  of  maximum 
length  in  its  direction. 

Conversely,  to  show  that  T)  A,  first  consider  the  case  where  K 
is  a  strictly  convex  body.  For  each  u  G  consider  the  points 

Xu  G  Hu  n  dK  and  X-u  G  H^u  C  dK.  The  continuous  function 
/  :  ^  defined  by 


fiu) 


X  'll  X _ 7, 


\X..  —  X- 


u\\2 
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has  the  property  that  {f{u),u)  >  0  for  all  u.  In  other  words,  the 
angle  between  u  and  f{u)  is  less  than  But  it  is  an  elementary 
exercise  in  topology  (see,  for  example,  Munkres  |12].  problem  10, 
page  367)  that  if  /  :  — >  5”*“^  maps  no  point  x  to  its  antipode 

—X,  then  /  has  degree  1  and,  in  particular,  is  surjective.  To  show 
that  T>  C  let  (p,  q)  G  Dy  for  some  v  G  S'™  By  the  surjectivity 
of  /  there  is  u  G  5™“^  such  that  f{u)  =  v.  According  to  the  last 
sentence  of  the  previous  paragraph,  XuX-u  is  the  unique  longest 
chord  parallel  to  v.  Therefore  p  =  Xu  and  q  =  X-u  and  consequently 
(p,  q)  eAu- 

The  case  where  K  is  not  strictly  convex  is  treated  by  a  standard 
limiting  argument.  Given  a  vector  v  G  and  a  longest  chord 

pq  parallel  to  v,  we  must  prove  that  (p,  q)  G  A.  Since  K  is  the 
intersection  of  all  strictly  convex  bodies  containing  K,  there  is  a 
sequence  {K^}  of  strictly  convex  bodies  containing  K  with  the 
following  two  properties. 

1.  There  is  a  longest  chord  pkqk  of  Kk  parallel  to  u  such  that 
limfc^oo  \\Pk  -  qk\\2  =  \\P  -  Qh  >  the  limits  lim^^ooPfe  =  p  €  K 
and  limfc^oo  Qk  =  Q  ^  K  exist. 

By  the  result  for  the  strictly  convex  case,  there  is  a  sequence 
of  vectors  Uk  G  5'™“^  such  that  pk  =  Kk  H  Hu^{Kk)  and  qk  = 
Kk  n  H-Uk{Kk)-  By  perhaps  going  to  a  subsequence 

2.  lim/j^oo  Uk  =  u  &  gm-i  gxists. 

It  follows  from  item  1  that  ||p  — g||2  =  ||7'“9|l2  p  —  q  is 
parallel  to  v.  Therefore,  pq  as  well  as  pq,  are  longest  chords  of  K 
parallel  to  f.  It  follows  from  2  that  if  H  and  H'  are  the  hyperplanes 
orthogonal  to  u  through  p  and  q  respectively,  then  H  and  H'  are 
parallel  supporting  hyperplanes  of  K.  Therefore,  necessarily  p  ^  H 
and  q  G  H' ,  and  consequently  (p,  q)  G  Au  C  A.  □ 

Definition  6.5  (Non- Antipodal  ITS).  If  A  C  M™  is  a  convex 
body,  then  /  :  M™  — >  M™  is  non- antipodal  with  respect  to  K  if 
f{K)  C  K  and  {x,y)  G  A{K)  implies  {f{x),f{y))  ^  A{K).  If 
T  =  {M™; /i, /2, ..., /at}  is  an  iterated  function  system  with  the 
property  that  each  /„  is  non-antipodal  with  respect  to  K,  then  T 
is  called  non-antipodal  with  respect  to  K. 

The  next  proposition  gives  the  implication  (4)  (5)  in  Theo¬ 

rem  [LTJ  The  proof  is  clear. 
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Proposition  6.6  (A  Topological  Contraction  is  Non- Antipodal).  If 
T  =  {K™’; /i, /2, /at}  is  an  affine  iterated  function  system  with 
the  property  that  there  exists  a  convex  body  K  C  such  that 
fn{K)  C  int(K)  for  all  n  =  1,2,  then  T  is  non-antipodal 

with  respect  to  K . 

The  next  theorem  provides  the  implication  that  (5)  ^  (1)  in 
Theorem  11.11 

Theorem  6.7.  If  the  affine  IFS  T  =  (K™'; /i, /2)  •••, /v)  is  non- 
antipodal  with  respect  to  a  convex  body  K ,  then  T  is  hyperbolic. 

Proof.  Assume  that  A'  is  a  convex  body  such  that  /  is  non-antipodal 
with  respect  to  K  for  all  f  £  IF.  Let  C  =  K  —  K  and  let  f{x)  = 
Lx  a  (z  IF,  where  L  is  the  linear  part  of  /.  By  Proposition  15.21 
the  set  C  is  a  centrally  symmetric  convex  body  and 

L(C)  =  L{K)  -  L{K)  =  f{K)  -  f{K)  <ZK-K  =  C. 

We  claim  that  L{C)  C  int{C).  Since  C  is  compact  and  L  is 
linear,  to  prove  the  claim  it  is  sufficient  to  show  that  L{x)  ^  dC 
for  all  X  E  dC.  By  way  of  contradiction,  assume  that  x  E  dC 
and  L[x)  E  dC.  Then  the  vector  x  is  a  longest  vector  in  C  in  its 
direction.  Since  x  £  C  =  K  —  K  there  are  xi,X2  E  dK  such  that 
X  =  X1—X2,  and  (xi,  X2)  E  F>{K)  =  A{K),  where  the  last  equality  is 
by  Theoreml6.4l  So  (xi,  X2)  is  an  antipodal  pair  with  respect  to  K. 
Likewise,  since  Lx  is  a  longest  vector  in  C  in  its  direction,  there  are 
2/1, 2/2  S  dK  such  that  Lx  =  yi  -  2/2,  and  (2/1, 2/2)  G  Td{K)  =  A{K). 
Therefore 

f{x2)  -  /(xi)  =  L(x2)  -  L(xi)  =  L(x2  -  xi)  =  Lx  =  yi-  2/2, 

which  implies  that  (fnixi),  fn{x2))  £  Td{K)  =  A{K),  contradicting 
that  /  is  non-antipodal  with  respect  to  K. 

If  dc  denotes  the  Minkowski  metric  with  respect  to  the  centrally 
symmetric  convex  body  C,  then  by  Proposition  15.41  C  is  the  unit 
ball  centered  at  the  origin  with  respect  to  this  metric.  Since  C  is 
compact,  the  containment  L[C)  C  int{C)  implies  that  there  is  an 
a  E  [0, 1)  such  that  IlLxll,^.  <  a  \\x\\(j  for  all  x  E  M"^.  Then 

dc{f{x),  f{y))  =  \\f{x)  -  f{y)\\c  =  \\Lx  -  LyW^ 

=  \\L{x  -  y)\\c  <  a  ||x  -  2/||c-  =  adc{x,y). 
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Therefore  dc  is  a  metric  for  which  each  function  in  the  IFS  is  a 
contraction.  By  Proposition  15.61  dc  is  Lipschitz  equivalent  to  the 
standard  metric.  □ 

7.  An  Answer  to  the  Question  of  Kameyama 

We  now  turn  to  the  proof  of  Theorem  11.21  the  theorem  that  set¬ 
tles  the  question  of  Kameyama.  If  A  C  and  T  =  (A;  /i,  /2, ...,  /at) 
is  an  IFS  on  A,  then  the  definitions  of  coding  map  and  point-fibred 
for  T  are  exactly  the  same  as  Definitions  12.41  and  12.51  with 
replaced  by  A.  The  proof  of  Theorem  11.21  requires  the  following 
proposition. 

Proposition  7.1.  //AC  and  T  =  (A; /i, /2, ..., /at)  is  an 
IFS  with  a  coding  map  vr  :  D  — such  that  7r(D)  =  A,  then  IF  is 
point-fibred  on  A. 

Proof.  By  Definition  12.51  we  must  show  that  limfc^oo  fai  °  fa2  ° 

■  ■■  o  ffj^{x)  exists,  is  independent  of  a:  G  A,  and  is  continuous 
as  a  function  of  cr  =  criCJ2  •  •  •  G  D.  We  will  actually  show  that 
limfc^oo  /(71  o  o  •  •  •  o  fauix)  =  vr(cr). 

Since  vr  is  a  coding  map,  we  know  by  Definition  12.41  that  fn  o 
7r((T)  =  TToSn{(y),  for  all  re  =  1,  2, . . . ,  A.  By  assumption,  if  x  is  any 
point  in  A,  then  there  is  a  r  G  D  such  that  vr(T)  =  x.  Thus 

lim  f^^o  ■  ■  ■  o  fau{x)=  lim  /afc(vr(r))  (since  7r(r)  =  x) 

k^oo  /c— >oo 

=  lim  7r(S(ji  o  8^2  o  ...  o  So-^.  o  t)  (by  Diagram  I2.1|] 

/c— >oo 

=  7r(  lim  So-j  o  3^2  o  ...  o  o  r)  (since  vr  is  continuous) 
k^oo 

=  7r(cj). 

□ 

Theorem  7.2.  If  F  =  (M™'; /i, /2, ..., /n)  is  an  affine  IFS  with  a 
coding  map  tt  :  D  — >  A,  then  F  is  point-fibred  when  restricted  to 
the  affine  hull  of  7r{Q).  In  particular,  if  ti{Q)  contains  a  non-empty 
open  subset  o/M™,  then  F  is  point-fibred  on  M™. 

Proof.  Let  A  :=  7r(D).  Since  /n(A)  C  A  for  all  re,  the  restriction  of 
the  IFS  F  to  A,  namely  F\a  '■=  {A;  /i,  /2, . . . ,  /at),  is  well  defined. 

It  follows  from  Propositionl/.llthat  F\a  is  point-fibred  and,  because 
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the  coding  map  for  a  point-fibred  IFS  is  unique, 


7r(cr)  =  lim  o  o  •  •  •  o  /<Tfe(«) 

/c— >oo 

for  (cj,  a)  &  Q  X  A.  It  only  remains  to  show  that  the  restriction 
•^lafr(A)  •=  (aff(^);  fi,  f2,  ■  ■  ■ ,  In)  of  the  affine  IFS  !F  to  the  affine 
hull  of  A  is  point-fibred. 

Let  X  G  aff(^),  the  affine  hull  of  A.  It  is  well  known  that  any 
point  in  the  affine  hull  can  be  expressed  as  a  sum,  x  =  for 

some  Ao,  Ai, Am  G  K  such  that  Yl^=o  uq,  oi,  G 

A.  Hence,  for  (cr,  x)  G  H  x  aff(H), 


lim  o  o 

/c— >oo 


lim  o  o  •  •  •  °  f<7kiY] 

/c— >oo 

p=i) 

m 

lim  V  Xpf^^  o  ■  ■  ■  o  (op) 

/c— ^•oo 

p=0 

m 

^  Ap7r(u)  =  7r(cj). 

p=0 

□ 


Theorem  11.21  now  follows  easily  from  Theorem  17.21  and  Theo¬ 
rem  11.11 

Proof,  fof  Theoremll.2p  Let  A  :=  7r(H)  and  let  dim  aff(H)  =  k  <  m. 
It  is  easy  to  check  from  the  commuting  diagram  12. II that  f{A)  C  A 
for  each  /  G  JF  implies  that  /(aff(H))  C  aff(H)  for  each  f  G  P. 
Since  aff(H)  is  isomorphic  to  Theorem  11.11  can  be  applied  to 
the  IFS  J^\aS{A)  •=  /ii  f2,  ■■■fN)  to  conclude  that,  since  it  is 

point-hbred,  is  also  hyperbolic.  □ 

Note  that  the  IFS  (M;  /),  where  f{x)  =  2x  +  1,  is  not  hyperbolic 
on  M,  but  it  is  hyperbolic  on  the  affine  subspace  {—1}  C  M. 


8.  Concluding  Remarks 

Recently  it  has  come  to  our  attention  that  another  condition, 
equivalent  to  conditions  (1)  —  (5)  in  our  main  result,  Theorem  1.1, 
is  (6)  T  has  joint  spectral  radius  less  than  one.  (We  dehne  the 
joint  spectral  radius  (JSR)  of  an  affine  IFS  to  be  the  joint  spectral 
radius  of  the  set  of  linear  factors  of  its  maps.)  This  information  is 
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important  because  it  connects  our  approach  to  the  rapidly  growing 
literature  about  JSR,  see  for  example  0,  0,  and  works  that  refer 
to  these. 

Since  Example  3.3  and  the  results  presented  by  Blondel,  Theys, 
and  Vladimirov  m  indicate  there  is  no  general  fast  algorithm 
which  will  determine  whether  or  not  the  joint  spectral  radius  of  an 
IFS  is  less  than  one,  we  feel  that  Theorem  1.1  is  important  because 
it  provides  an  easily  testable  condition  that  an  IFS  has  a  unique 
attractor.  In  particular,  the  topologically  contractive  and  non- 
antipodal  conditions  (conditions  4  and  5)  provide  geometric/ visual 
tests,  which  can  easily  be  checked  for  any  affine  IFS.  In  addition 
to  yielding  the  existence  of  an  attractor,  these  two  conditions  also 
provide  information  concerning  the  location  of  the  attractor.  (For 
example,  the  attractor  is  a  subset  of  a  particular  convex  body.) 
We  also  anticipate  that  Theorem  1.1  can  be  generalized  into  other 
broader  classes  of  functions,  where  the  techniques  developed  for  the 
theory  of  joint  spectal  radius  will  not  apply. 
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